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Abstract 

The review deals with the physics of cold atomic gases in the pres- 
ence of disorder. The emphasis is on the theoretical developments, 
although several experiments are also briefly discussed. The review is 
intended to be pedagogical, explaining the basics and, for some of the 
topics, presenting rather detailed calculations . 

1 Introduction 

The behavior of cold atoms and Bose-Einstein condensates (BEC) in the 
presence of a random potential is an active area of research, and several re- 
views on the subject had already appeared [H [21 El IH El E]- Those reviews 
discuss, among other things, the experimental work on BEC spreading in 
a disordered wave guide, which have culminated in the observation of one- 
dimensional Anderson localization [71 18]. The present review is devoted to 
two and three dimensional systems of cold atoms, and it concentrates on top- 
ics not extensively covered in the previous reviews. There exists considerable 
amount of theoretical work on the subject, and recently few experimental 
results on spreading of an atomic cloud, in two [9j and three [TOl [TT] dimen- 
sions, have been published. The review is intended to be self contained and 
pedagogical, accessible to people without a solid background in the physics of 
disorder. It therefore contains a significant amount of background material. 

Disorder is, of course, a very old theme in condensed matter physics, and 
a lot of knowledge was accumulated over the years. Some of this existing 
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knowledge can be directly transferred to the field of cold atoms, although 
there are also considerable differences between the two fields. For one thing, 
a random potential for atoms is usually created optically (an optical speckle) 
and it differs, in some respects, from the common models used in condensed 
matter theory. Furthermore, in disordered electronic systems one is usually 
interested in the nature of eigenstates at the Fermi level, because this de- 
termines whether a given system is a metal or an insulator. In contrast, 
a typical transport experiment with an atomic BEC involves an expansion 
of the atomic cloud ( a "matter wave packet") in the presence of disorder. 
Such a packet contains a broad spectrum of components, some of which can 
propagate whereas the others get locahzed. Nonhnear effects also come into 
play, especially at an early stage of the expansion. 

The equilibrium properties of a disordered BEC are also of interest, in 
particular the phenomenon of the insulator-superfluid transition. When the 
density of particles is low, they can at most form a fragmented condensate, 
i.e. small disconnected lakes, with no coherent coupling between them (an 
insulator). With an increase of the particle density these lakes overlap and, 
at some critical density, a single coherent condensate is formed (a superfluid). 

It would hardly be possible to cover, in a single review, all aspects of the 
broad and rapidly progressing field of disordered cold atom systems. The 
choice of topics reflects, to some extent, personal preferences of the author, 
and every effort was made to cover those topics with some depth, clarity and 
precision. The emphasis is made on identifying the relevant parameters and 
on discussing limiting cases, when simple and clear physics emerges. Such 
simplified discussion might not be sufficient for quantitative interpretation of 
the existing experiments but, at least, it provides a good starting point for 
appreciating the richness and complexity of the underlying physics. 

Sections 2 and 3 of this review contain background material on speckle 
potentials and on behavior of a quantum particle in a random potential. 
Section 4 deals with a disordered BEC in equilibrium and discusses the con- 
ditions for a disorder induced insulator-superfluid transition. Section 5 is 
devoted to the dynamical problem of cold atoms expanding in the presence 
of disorder. Both BEC and Fermi gas are treated. A dynamical problem of a 
different type is considered in Section 6. There a BEC with an initial phase 
modulation undergoes free expansion. In the course of the expansion large 
density fluctuations can develop, resulting in "matter wave caustics". Some 
concluding remarks are presented in Section 7. 
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2 Optical speckles: A random potential for atoms 



An external potential (not necessarily random) for atoms can be created by 
lasers |12t [T3]. The laser radiation polarizes an atom and changes its energy, 
e.g. the ground state energy, by an amount 

V^(r) = ^«MJ(r), (1) 

where /(r) is the radiation intensity at the position r of the atom and a{uj) 
is the atomic polarizability. The latter depends on the radiation frequency co 
and on the properties of the atom, in particular, on the detuning 6 = u — Ur, 
where ojr is a resonance frequency of the relevant atomic transition. The 
case 6 > (blue detuning) corresponds to a positive (with respect to the 
unperturbed atomic ground state) potential, while 6 < (red detuning) 
corresponds to a negative potential. Thus, a random pattern of light in- 
tensity, referred to as speckle pattern, will serve as a random potential for 
atoms. Such speckle potentials are widely used in experiments with cold 
atoms [71 ini Uni [m [in Uni Uni [TT] . Although there exists a rather compre- 
hensive book on optical speckles jTH], it is worthwhile to review the subject 
here, paying attention also to the most recent work [HI [201 El 122] • 

2.1 Qualitative discussion 

A schematic set up for making a speckle pattern is depicted in Fig.l. 
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Figure 1: A laser beam, passing through a diffuser, illuminates an area of size 
D at the outgoing plane z = (the source plane). Light from an elementary 
source of size w (three such sources are designated by black spots) propagates 
in a cone of angular width \/w and illuminates a region of size z\/w in 
the observation plane. Waves emerging from different elementary sources 
interfere with each other, creating a speckle pattern. 

A laser beam is transmitted through a diffuser. A common example of a 
diffuser is a glass plate with a rough surface (ground glass). Such a diffuser 
imprints a random phase on the transmitted radiation, while the intensity at 
the outgoing plane (the plane 2 = 0) simply follows the shape of the beam, 
e.g., a Gaussian laser beam. In the literature [18j one often studies a geometry 
when a plane wave (the central part of a laser beam) impinges on a circular or 
rectangular hole in a opaque screen, the hole being covered by a ground glass 
diffuser. In this case the light intensity is uniform within the corresponding 
circular or rectangular domain, and it is zero outside that domain. The 
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detailed geometry of the set up is not essential for the qualitative arguments. 

Thus, an illuminated area, of size D, appears at the outgoing plane of the 
difFuser (the z = plane, often referred to as source plane). This illuminated 
area can be viewed as a collection of sources, of a typical size w (usually, of 
the order of the radiation wavelength A) , with no spatial coherence between 
different sources. The sources are monochromatic, of frequency cu, and the 
n'th source has a phase factor e-^(4>n+i^t) ^ Here n runs from 1 up to the total 
number of sources, ~ (D/w)^, and the phases 0„ change randomly from 
one source to another. The collection of such random sources at the z = 
plane will produce a random, three-dimensional intensity pattern ("grains" 
of light) in the z > half-space. By choosing some distant observation 
plane, z = const, one can study the transverse properties of such speckle 
patterns. Their longitudinal properties, such as the typical grain size in the 
^-direction, are also of interest. Although the detailed theory of speckles is 
rather intricate, a rough qualitative picture can be obtained by using a few 
basic facts about interference and diffraction. 

We discuss the transverse properties first. An elementary source, of size 
Wj in the source plane can be thought of as a particle of this size which 
scatters the incoming light or, alternatively, as a pinhole which diffracts the 
light. In either case the light will propagate in a cone of angular width 
~ X/w, illuminating a region of size ~ zX/w in the z-plane. Let's consider 
first the case {zX/w) > D, i.e. z > (Dw/X) = z*. This case corresponds to a 
"fully mixed" speckle, when all elementary sources contribute to the intensity 
at a typical point P in the 2;-plane. The electric field at some point P in 
the observation plane is a sum of many complex amplitudes, with random 
phases, originating from different areas at the source plane. As a result, the 
light intensity / across the observation plane is a random function of the 
position P, with the characteristic bright and dark spots. An estimate of the 
typical size of such speckle spot is analogous to the estimate of the separation 
between neighboring maxima and minima in the standard diffraction theory 
(the difference is that in the diffraction pattern from, say, a circular hole 
maxima and minima occur only in the radial direction, whereas in the speckle 
pattern they will be encountered in all directions). Let us denote by ap the 
phase difference between the two waves arriving to the point P from the 
opposite edges of the illuminated area at the source plane. If the point P is 
moved by a distance s±, the phase difference ap will change by ~ s_iD/Xz, 
i.e., a shift of the observation point by a distance Sj_ ~ Xz/D will change 
constructive interference to destructive. The distance s±_ provides an estimate 
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for the size of one speckle spot. Note that our discussion pertains to what is 
called "a free space geometry" [18j, when the light propagates freely between 
the source plane and the observation plane. In practice one usually employs 
the "imaging geometry" [18], when a converging lens is placed between the 
two planes and the observation plane coincides with the focal plane of the 
lens. The speckle spot size then becomes s± ~ Xf/D, where / is the focal 
length of the lens, and the minimal size is determined by the diffraction limit. 

The picture outlined in the previous paragraph undergoes a profound 
change when the distance z to the observation plane becomes smaller than 
z* - the so-called deep Fresnel zone [121 [201 EH 122] • According to the above 
estimate, for z ^ z* a speckle spot reaches a size s± ~ w. When z decreases 
further, a point P in the 2-plane can receive radiation only from a limited 
number of the elementary sources, namely, the size of the region in the source 
plane which contributes to illumination of point P is -De// ~ Xz/w. Using 
the above estimate for the speckle spot size, with -De// instead of D, one 
obtains s± ~ Xz/D^^ff ~ w. Thus, s± ^ w holds starting from z ^ z* and 
all the way down to the minimal distance z ~ w'^/X. 

We now briefly discuss the structure of a speckle pattern in the longitu- 
dinal direction. Also in this direction the light intensity exhibits an irregular 
sequence of maxima and minima, and separation between those defines the 
typical size, s\\, of the light "grain" in the z-direction. Standard diffraction 
theory [23] tells us that a significant change of intensity in the 2;-direction 
will occur when the corresponding phase, -D^/ Xz, will change by amount of 
order unity, i.e. such a change will occur over a distance s\\ ~ A(z/-D)^. In 
the deep Fresnel zone, z < z*, one has to substitute D^fj for D, obtaining 

[211 122] s|| ~ wVA. 

2.2 Analytic treatment 

Neglecting possible polarization effects, we consider a scalar field which can 
be thought of as one particular component of a monochromatic electric field. 
We denote the field at the source plane hy ijj{p,0), where the vector p desig- 
nates the in-plane position and the second argument stands for z = 0. The 
function tp{p,0) is significant only within the region \p\ ~ D. The position 
in the observation plane z is denoted by R. The field in this plane created 
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by the source field ip{p, 0), is given by 



(27r)2 

where = u / c = 271 / \, and 



Vi(q,0)= / dpij{p,^)e-'^-P. 



(2) 



(3) 



The expression Q for V^(R, 2;) is exact, containing both propagating and 
evanescent waves, and it has a simple interpretation. For z = Q the source is 
Fourier expanded with respect to the transverse coordinate p, and the factor 
exp(iA/A;o ~ l"^^) "Fourier propagator" which relates the transverse 

Fourier transform at the z-plane, ?/'(q, z), to '?/'(q, 0). 

In the paraxial approximation the Fourier propagator is approximated 
by e-x.Yi[ikQZ — {iq^z/2kQ)]. Substituting ([s]) into ([2]), with the approximated 
Fourier propagator, and performing integration over q leads to: 



^(R,z) 



'^0 ^ikpz 



2Triz 



dpipip, 0) exp 



•^0 ^-D 



(4) 



which is the mapping of the source field, ilj{p,0), onto the field '?/'(R, 2;) at 
the observation plane. 

ijj{p, 0) is a random function of the position p. It is described by the 
correlation function tp{p, 0)tp*{p', 0), where the overbar indicates a statistical 
average. Eq. Q enables one to immediately write down the correlation 
function for the bulk speckle pattern, tp{R, z)tp*{R', z'). We give the explicit 
expression only for the "transverse" correlation function {z = z'): 



i){R,z)i)*{R',z) 



kp 
2nz 



dpdp'ip{p, 0)iIj*{p', 0) exp 



{AR-Ap)-{Ro-po) 

Up + p')- 



(5) 



where AR = R- R', Ap = p- p', and Ro = l{R + R'), po 

This expression simplifies if one assumes a white noise source, i.e. ipiP^ 0)V'*(p') 0) 
Xl{p)S{p — p'), where x is a proportionality constant and I{p) = \ip{p,0)\'^ 
is the average intensity at point p in the source plane. For this case one 
obtains 1181: 



^{R,z)ip*{R',z) 



X 



kn 



2'KZ 



exp 



t^Ro-AR 



dpl{p) exp 



-t^AR-p 

z 



(6) 
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which is proportional to the Fourier transform of I{p). If /(p) signifi- 
cantly differs from zero in a domain of size then the correlation function 
iplR, z)iIj*{R' , z) decays on a characteristic scale s± ~ Xz/D, which defines 
the typical transverse size of a speckle spot, as estimated in Sec. 2.1. For 
simple cases, when I{p) is constant within an aperture of a rectangular or 
circular shape, the integral in ^ can be evaluated explicitly P^8J. Another 
simple and instructive example is I{p) = Jq exp(— p^/2D^). The integral in 

(6) is then given by 27rD^exp(— ^ 22^^^ ) 'which, again, reveals the charac- 
teristic decay length Xz/D. 

One can relax the assumption of a white noise source and consider a more 
general correlation function, at the source plane 



^{p,0)r{p',0) = C{Ap)I{po), (7) 

where the function C decays on a shot scale, w, whereas the average intensity 
/ changes on a much larger scale, D. Using this expression as an input in 
Eq. ([5]), one can derive the transverse properties of the speckle pattern. The 
term kQ(p ■ Ap)/z ~ koDw/z, in the exponent of Eq. ([s]), can be neglected 
for z > Dw/X = z*. In this case the double integral in ([s]) factorizes into 
a product of two integrals which are the Fourier transforms of /(po) and 
C{Ap). These Fourier transforms determine, respectively, the speckle spot 
size s_L ~ Xz/D and the extent Xz/w of the entire speckle pattern in the 
2;-plane. In the opposite case, z < z*, the speckle pattern "freezes", as 
described in Section 2.1. 

The more general correlation function, ip{R, z)iIj*{R', z'), with z 7^ z', has 
been studied in detail in |21||22]. where the normalized function, 7(-R, z; R', z') 



ip{R, z)ip*{R', z') / y I{R, z)I{R', z') was introduced (this function is denoted 
there by /i). Simple analytic expressions for |7(-R, z; R + AR, z + A^;)^ can 
be obtained if one chooses in ([7| C{Ap) = exp[— (Ap/w)^] and /(po) = 
Jo exp[— (po/-D)^] [211 [22]- In certain limiting cases the expression 



\j{R, z;R+AR,z + Az)f = ^ exp 

1 + (A2;/S||)^ 



2{AR/s^f 



l + (Az/s||)^ 

emerges. In the deep Fresnel zone = nw'^/X and s± = w (Eq. (12) in |22]). 
whereas in the Fresnel zone (and for Az not too large) sy = Xz'^/nD'^ and 
s± = Xz/nD |25| . The correlation function (|8| is quite unusual: not only is it 
strongly anisotropic but also the effective correlation length in the transverse 
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direction, s±eff ~ s±^,/Y^^~(AzJs^ , depends on the separation Az between 
the points in the longitudinal direction. Moreover, s±eff increases with Az, 
so that long range correlations in the transverse direction develop. 

A simpler three-dimensional potential, corresponding to the interference 
pattern in an ergodic cavity, was considered in It is fully isotropic and 
the square of the correlation function is 

^^(r) = smc^(27rr/A), (9) 

where A is the light wave length. 

Let us recall that the potential acting on the atoms is proportional not 
to the field but to its intensity, see Eq.Q. Therefore, denoting {R, z) = 
r, the correlation function for the potential, V{r)V{r'), is proportional to 
I{r)I{r') = %l}{r)ip*{r)'il){r')ip*{r'), which requires averaging of a product of 
four fields. For a Gaussian random function the latter reduces to the sum of 
two possible contractions of the fields, with the result |18] : 



I{r)I{r') = I{r)I{r') + \^{r)ilj*{r')\^. (10) 

Since /(r) changes on a scale much larger than the decay length of the 
correlation function, one can set /(r) = I{r') = Jq. Rewriting (10) in terms 
of the potential then yields j^Hl 



V{r)V{r') = V,'[l + \^{r-r')\% (11) 



where Vq is the average value of the potential at point r and 7 = il){r)i})*{r') / 1^ 



is the normalized correlation function of the speckle field. It follows from ( 11 ) 



that the second moment = 2V^. The n'th moment, V^, will involve av- 
eraging of n pairs of fields and thus n\ contractions, giving = n\V^. This 
implies that the probability distribution of the potential value, at some point 
r, is given by the Rayleigh distribution 

P(V) = ^expi-V/Vo), (12) 

where a blue detuned potential has been assumed (for red detuning both V 



and Vo are negative, so that (12) should be written with a minus sign). The 
Rayleigh distribution for the intensity is a direct consequence of the fact that 
the speckle field - being a sum of many contributions with random phases - 
obeys the central limit theorem, i.e., the real and the imaginary part of the 
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field are Gaussian distributed [IB]. The distribution PiV) is not symmetric 
with respect to its average value Vq, and it is convenient to introduce the 
deviation 6V = V — Vo, with the correlation function 



where I7P was denoted by F and the argument (r — r') was rescaled by 
the correlation decay length i?o- Computer generated images of a speckle 
potential (in 2d) can be found, e. g., in |26| 128]. 

Finally, let us mention that making an optical speckle is not the only way 
to create a random potential for atoms. Another option is to introduce into 
the system some foreign, "impurity" atoms, trapped in a random fashion at 
the nodes of an optical lattice [29j. These impurity atoms form a frozen, 
random pattern which scatters the mobile atoms of the system under investi- 
gation. This method for creating a random potential has been implemented 
in the recent study of a one-dimensional Bose system [30j. 

3 Quantum particle in a random potential 

The behavior of a quantum particle in a random potential has been exten- 
sively studied in condensed matter physics [211 1231 EH EH EHl E3- The 
theoretical tools range from diagrammatic technique [35], to self-consistent 
approach [38j, to field theory methods [36, 37J. The most popular model is 
that of a Gaussian random potential, with or without correlations (the latter 
case is known as the "white-noise limit"). This model features not only in 
condensed matter physics but also in the field of cold atoms, for instance, 
in the study of the disorder induced superfiuid-insulator transition [39| HO] . 
The Gaussian potential differs in some respects from a speckle potential. The 
latter is bounded from below (blue detuning) or from above (red detuning). 
This will clearly lead to large changes in the particle behavior close to the 
boundary of the spectrum. In addition, the speckle potentials lack symme- 
try (in statistical sense, of course), with respect to the average value, so that 
products of an odd number of 6V^s do not average to zero. Therefore the 
diagram technique for a speckle potential [261 127] differs somewhat from the 
standard one [35]. To a large extent, however, there are no qualitative dif- 
ferences, as long as one is in the weak scattering regime. In particular, the 
mean free path, calculated in the Born approximation, is "universal", in the 




(13) 
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sense that it only depends on the correlation function (13) but not on the 
detailed statistics of the potential. 

Consider a random potential V{r), with some average value V. This 
value is included into the unperturbed Hamiltonian, and the perturbation 
is SV{r). Thus, the bottom of the unperturbed spectrum is V and, if the 
particle energy e is reckoned from this point, the unperturbed spectrum is 
simply e = Jr^k"^ /2m. With this caveat, most of the ensuing discussion applies 
to a broad class of potentials, including speckle and Gaussian potentials. 
Viewing the potential relief from the "V^-level", one can envisage a collection 
of potentials barriers and wells, of typical magnitude Vq and typical size Rq. 
It is convenient to define the "correlation energy" 

where m is the particle's mass. Eq has the clear meaning of a zero point 
energy for a particle confined to a spatial region of size Rq. There are, 
thus, three relevant energy scales in our problem - the typical variation Vq 
of the random potential, the correlation energy and the particle energy 
e - and the physics depends on the two ratios: Vq/Eq and e/E^ = [kRo)"^. 
In particular, kRo << 1 corresponds to short range correlations, when the 
correlation length Rq is shorter than the de Broglie wavelength of the particle. 
The opposite case, kRo » 1, describes smooth (on the wavelength scale) 
disorder, for which a semiclassical approximation might be a good starting 
point. The space dimension, d, is also an important factor. For d = 1,2 all 
eigenstates (in an infinite system) are localized, while for d = 3 there is a 
critical energy Ec which separates between extended and localized states. In 
this review we concentrate on two and three dimensions. 

The value of the parameter (Vq/Eq) = rj is crucial for understanding the 
transport picture. For rj << 1 a typical well is too weak to bind a particle. 
In this regime semiclassical considerations are not appropriate and classical 
percolation is of no relevance whatsoever for understanding the quantum lo- 
calization transition. In the opposite case, rj » 1, there is a broad region 
in the particle energy spectrum, where semiclassical considerations do apply. 
The ultimate transition to localization is, of course, of quantum nature also 
for this case but now the transition occurs close (somewhat above) to the 
classical percolation threshold and an interesting crossover between classical 
percolation and Anderson localization takes place. All this is discussed below 
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in some detail. Before making various estimates and calculations, it is essen- 
tial to specify whether 7] is large or small, and it is convenient to separately 
discuss the following cases: 



3.1 d = 3 and 77 << 1 

A particle of sufficiently high energy will be weakly scattered by the random 
potential and will propagate by diffusion. When the energy of the particle 
is lowered, the scattering becomes stronger and quantum interference effects 
cease to be negligible. These effects lead to Anderson localization |41j below 
some energy Ec, called the mobility edge. 

Before trying to estimate Ec let us make the following remark: The 3d 
speckle potential is not "generic" in the sense that the integral j r{r)d^r 
diverges, due to the long-range character of the correlation function r(r) (see 
equations ([s]) and ([9])). In contrast, the impurity potential, mentioned at the 
end of Sec. 2.2, belongs to the standard type of a "generic" potential, usually 
encountered in condensed matter physics. The expression for the mean free 
path, and therefore an estimate for Ec, depends on whether the potential is 
"generic" or not. In the forthcoming discussion the "generic" case will be 
considered first and, then, some necessary changes for the speckle potential 
will be pointed out. 

Let us now estimate the characteristic energy near which the diffusive 
transport breaks down. Qualitatively, the random potential can be viewed as 
a collection of scatterers, of strength Vq and size Rq. Since Vo « Eq, a. single 
scatterer can be treated within the Born approximation, regardless of the 
particle energy. One should, however, distinguish between "slow" particles, 
with energies e << Eq (i.e. fci?o << 1) ^"^^ "fast" particles, satisfying 
e >> ii^o (i-e. kRo » 1). For slow particles the scattering cross-section is 
estimated as |42) a ~ rj^R^- Since the number of scatterers per unit volume 
is Uscat ~ Rq^ 1 the mean free path 



When one refers to "a quantum particle" , one actually has in mind a wave 
packet with a well defined wave vector k, such that the uncertainty A/c < k. 
The minimal size of such a packet is of the order of k~^ and, for the above 
reasoning to be consistent, this must be much smaller than /, which leads to 




(15) 
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the requirement kl » 1. The condition kl ^ 1 (the loffe-Regel criterion) 
corresponds to a characteristic energy 

W-V,v'-{^y{V,^Rl)\ (16) 

Thus, when the particle energy e (counted from the average value of the 
potential) decreases and approaches W, the diffusive transport breaks down 
and quantum interference effects come into play in an essential way . The 
potential strength and range appear in Eqs. ( T5p6 ) only in the combination 

We now return to the fast particles, with e >> Eq. Such high energy par- 
ticles scatter primarily in the forward direction, within a cone of an angular 
width < l/kRo- Therefore the scattering cross section for fast particles 
is (7 ~ rf/k'^ |32j (smaller by a factor [l/kR^Y than the cross section for 
slow particles). Furthermore, the cross section relevant for transport is not 
a but the transport cross section, a* ~ a/lkRo)"^. The small extra factor, 
[l/kRoY, accounts for the fact that in a single scattering event the particle 
direction changes only by a small angle, randomly distributed in a cone of 
angular width A6 ~ l/kR^. Therefore it will take [kRo)"^ scattering events 
to completely randomize the initial direction. Thus, the transport mean free 
path is: 



The qualitative estimate in Eqs. (15), (17) can be supported by a calcu- 



lation based on a diagrammatic perturbation theory [35j . In that theory the 
transport mean free time r* is related to the correlation function of the po- 
tential. In the leading order in Vq, the inverse mean free time (the scattering 
rate) is expressed in terms of the Fourier transform r(q) of the correlation 



cos ( 



function r(;^), defined in Eq. (13): 

T* h ^ J (27r)3 ^ ^ V 2m 

= Jl^o'' desinet (^2ksm^^ (1 -cos^), (18) 

where e = is the particle energy and u = mk is the density 
of states (per spin) at that energy. The argument 2k sin | corresponds to 
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the momentum transfer for (elastic) scattering at an angle 9 and the factor 
(1— cos 9) accounts for the difference between the scattering and the transport 
cross sections. Taking, for example, a Gaussian correlation function, F = 
exp(— i?^/i?Q), performing the integral and writing the result in terms of 
/* = VT* { V = hk/m is the particle velocity), yields 

2Rn ^ eV 1 



^* ^[voj i-ii + k^Dexpi-kmiY ^^^^ 

For kRo << 1 one can expand the exponent, obtaining I* = ARQ/rfy^ 



(compare to the qualitative estimate in (15)), whereas in the opposite case 
the result is simply /* = ^^(^)^ (compare with (17)). Eq.(19) can be re- 
written in terms of the diffusion coefficient = vl* 
the result is: 



/S^ For thehmiting cases 



-T^fn^l , W «{e-V) «Eq 



1/2 

, W «[e 

— \ 5/2 

) ' ^ » ^0 

Note that we have restored V in these expressions, to give them the more 
convenient form, "invariant" with respect to the choice of the reference en- 
ergy. Other types of the correlation function F(R) can be considered in a 



similar way. The order-of magnitude estimates (15), (17) do not depend on 
the specific shape of the correlation function, provided that J F(R)(i^i? is 
finite. 

Thus, for the case r] « 1 the Born approximation is valid in a broad 
range of energies, the only condition being (e — V) » W . When (e — V) 
approaches W , the Born approximation breaks down and the diffusion coef- 
ficient becomes of the order of h/m. This is the natural unit of "quantum 
diffusion" (analogous to the quantum conductance e^/K) and it signals a 
crossover to a purely quantum mode of transport. At the diffusion coeffi- 
cient drops to zero. Somewhat above it obeys a power law, oc {e — EcY, 
with u ^ 1.57 Below Ec all states are localized. The exact position of 
Ec depends on the type of the random potential. It should be located some- 
where in the vicinity of the average potential value V, although the precise 
location, and even the sign (with respect to V), is not known. The aforemen- 
tioned loffe-Regel criterion marks breakdown of the diffusive transport but it 
does not enable one to determine the location of E^. In technical terms, the 
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mean free path is related to the imaginary part of the self-energy, whereas 
Ec can be affected also by the real part, which causes a shift of the spec- 
trum. The real part, calculated in the self-consistent Born approximation 
(see ^44], [45] and Sec. 5.3), turns out to be negative, of the order of rjVo, 
which indicates that Ec is below V. One should keep in mind, however, that 
the self-consistent Born approximation is not really a controlled one. 

For energies well below Ec, when {Ec — e) >> W, one enters the region 
of strongly localized states- the "Lifshitz tail". While for high energies the 
density of states z/(e) is close to that of a free particle, in the Lifshitz tail 
z/(e) is small and it depends on the type of the potential. For a Gaussian 
potential the tail was studied long ago, by the method of optimal fluctuation 
|16],I171[32]. At extremely low energies, [Ec — e] > Eq, the optimal fluctuation 
is just a single well, of size Rq and depth ~ |e|, which immediately leads to 
lnz/(e) oc —e'^/VQ. For energies W << Ec — e << Eq the optimal fluctuation 
is of different nature. In this energy interval a single well of size Rq and depth 
~ |e| cannot bind a particle, and the optimal fluctuation corresponds to a 
much broader well, of radius R^ ~ {h^ /m\e\yf'^ » Rq. Since such a well is 
comprised of ~ {R^/RqY "elementary" wells, the corresponding probability 
is exp(— c^/Vq^), raised to power {R^/RqY. This yields 



e 



2 \ / f^2 \d/2 ^, |x2-(d/2) .^4^l/(4-d) 



In 3d, Et = W, and in 2d it is the characteristic crossover energy Vq/Eq, 
which will be identified in Sec. 3.3. It is interesting to note that the same en- 
ergy scale W (in 3d) emerges both from the loffe-Regel criterion and from the 
analysis of the Lifshitz tail. The factor ~ {R^/RqY, entering the derivation 



of (21 ), can be also interpreted in a different way: Since the spatial extent of 
the localized wave function, R^, strongly exceeds the typical variation length 
Rq of the random potential, the particle experiences not the original poten- 
tial but an effective Gaussian potential with a variance smaller by a factor 

{rjrqY 



In 3d, (21) can be written as lnz/(e) ~ —C/R^, where C ~ Ro/v is a 
characteristic length, specifying the short range disorder, and termed the 
Larkin length in |l0] (it is of the same order of magnitude as the mean free 



path in the weak scattering regime, (15)). The localized states in the Lifshitz 



tail are rare, separated by exponentially large distances, and the tunneling 
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amplitude between a pair of such states is exponentially small. When e 
approaches the mobility edge, by a few Et, the size of the localized wave 
function approaches C, the overlap becomes significant and the tunneling 
amplitude is of the order of unity: the concept of the optimal fluctuation 
ceases to be relevant, 

So far there have been no quantitative studies of the Lifshitz tails for 
speckle potentials in 2d and 3c? (the Id case was considered in |1H])- For a 
red-detuned potential, which is not bounded from below, one can argue that 



Eq. (21) for the Lifshitz tail is still valid. Although the original potential is 



not Gaussian (see (12)), the particle localized in a optimal fluctuation feels 
an effective potential, obtained by integrating the original, "microscopic" 
potential over a volume Rf. This results in an effective potential with an 
approximately Gaussian statistics, with a variance by the factor {R^/R 







smaller than that of the original potential, so that (21) is recovered. This 
consideration does not hold in the extreme tail, |e| < Eq, when the optimal 
fluctuation is just a single well of depth ~ |e| and size ~ Rq. The density of 
states then follows the tail of the potential distribution (12), i.e. lni/(e) oc 
-\e\/Vo. 

For the blue-detuned speckle potential the picture is different. Here one 



can expect the validity of Eq. (21 ) only if the particle energy satisfles W « 
(Ec — e) << Vq, i.e. when e is well below E^ but still far away from the bottom 
of the spectrum. When e approaches the bottom of the spectrum, e = (here 
it is natural to reckon the energy from the bottom of the spectrum), a new 
tail, of different nature emerges. Close to e = the density of states can be 
estimated as follows [32] (the argument applies in any dimension d). In order 
to arrange for a state with a small energy, one has to "clean out" a region of 
size R^ ~ {h'^/me), i.e. no barriers of height larger than e should be allowed 
in that region. Such a region consists of many, ~ {R^/RqY, independent 
elements, of size Rq each. The probability that the potential in a given 
element does not exceed e is ~ e/Vo. Therefore, for e — )• 

v{e)(x{y\ ^ lnz/(e) oc - f In-j (22) 



The two expressions, (21 ) and (22 ) should match somewhere half way between 
the bottom of the spectrum and E^.. 

As has been stated in the beginning of this section, for the 3d speckle 



potential the integral jT{Y)d^r diverges, so that the expression (15) for the 
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mean free path and the subsequent estimate for W do not apply. The formula 
(18), however, does apply and, for the ergodic cavity potential ([9]), it yields 
for the mean free path (in the small k limit, when there is no difference 
between I and /*) 1261: 



k 



(23) 



with i?o = A/27r. Thus, for kRo << 1, the mean free path is proportional 
to k, unlike the fc-independent result in (15). The condition /c/ ~ 1 will now 



give a characteristic energy Vq/Eq, instead of (16). This energy provides an 
estimate for E^, with respect to V. 

Finally, the study of transport in a 3d anisotropic speckle potential, with 
correlation function ([s]) (or, possibly, even some more general and compli- 
cated cases |21l 122] ) only begins and the first work on the subject has just 
appeared jl9]. Due to the anisotropy, diffusion is described by a tensor which, 
in the weak scattering regime, can be computed with the help of the Fourier 
transform of 



r(q±,gii 



TTV 27r 



s±s\\ 

(1± 



exp 



qlsl 



q\\s\\ 

q±s± 



(24) 



where q_|_ and are the transverse and the longitudinal (i.e. z) component 
of the transmitted momentum. Because of the long range transverse correla- 
tions in the speckle potential (see the discussion after Eq. (8)) r(qj^, gy) has 
a peculiar essential singularity at = 0. More theoretical work is needed 
in order to understand how such strongly anisotropic, "non-generic" poten- 
tials affect the standard picture of Anderson localization. It is worthwhile 
to note, however, that the subtle long-range correlations, of the kind given 
in (|8]), although important in principle, might be not so prominent in prac- 
tice. For instance, according to the authors of [10|, their speckle potential is 
well described by a (anisotropic) Gaussian correlation function which is, by 
definition, "generic" and short-ranged. 



3.2 d = 3 and r/ >> 1 

This is the case of a "smooth" random potential (the condition Vq » Eq will 
be always satisfied if the correlation length Rq is made large enough). Since, 
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in contrast to the previous case, the energy Eq is much below the height of a 
typical barrier, the condition kRo » 1 becomes now the necessary condition 
for diffusive transport. We start with the Born approximation, which should 
always hold at sufficiently high energies. 

For the case Vq >> Eq, the Born approximation for a potential barrier of 
strength Vq and size Rq results in a scattering cross section 

"-(^y. (25) 



and the approximation is valid only if the particle energy satisfies |42] 

e»^= Ea. (26) 
Eq 

The scattering is strongly anisotropic and the transport cross section a* is 
smaller than a by a factor {l/kRo)"^, which yields a transport mean free 



path given in Eq. (17). The latter contains only quantities with a clear 
classical meaning, and it can be arrived at by purely classical considerations. 
Consider a particle with an initial velocity v, in a given direction. The 
particle is subjected to a force of magnitude |F| = |V^(r| ~ Vq/Rq- The 
direction of the force is changing randomly, at a distance of order Rq. During 
a time interval At = Rq/v the particle will acquire a velocity increment 
\Av\ ~ EAt/m ~ Vq/tjiv. This increment has a random direction, so that 
the number of time intervals At, needed to completely degrade the initial 
direction, is of order (f/Af)^ ~ (e/Vo)^. Thus, the time after which the 
initial direction is forgotten (the transport time) is 



r ~ 



f \ » / e \ R() 



and the corresponding mean free path is 



2 

e 



I* ^vr* ^Ro{-] , (e » Vo) (28) 



which coincides with (17). This coincidence was pointed out in [50j and it 



is quite remarkable, because Born approximation is incompatible with semi- 



classics [42] . Indeed, the Born scattering cross section, Eq. (25), strongly 



differs from the semiclassical one, which is of order R^. Also the differential 
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cross sections are entirely different: the typical scattering angle in the Born 
approximation is l/kRo, whereas in semiclassics it is Av/v ~ Vo/e. It is 
only the transport scattering cross section (and hence the transport mean 
free path) that is given by a similar expression in both approaches. Further- 
more, for e < the Born approximation breaks down. The estimate for 



/* in Eq. (28), however, holds as long as e >> Vq. Thus, we arrive at the 



conclusion that the only condition required for a diffusive transport (with 



(28) serving as an order-of-magnitude estimate for /*) is e >> Vq, i.e. the 
particle energy must well exceed the typical hight of the potential barriers. 
The corresponding diffusion coefficient is 

A = vr/3 ~ Doie/Vof/^ (29) 

where Dq = (Vq-Rq/'^)^^^ the natural unit for classical diffusion (this is 
the only combination, of the three quantities involved, with dimensions of a 
diffusion coefficient). 

When e approaches Vq, the transport mean free path approaches a value of 
order Rq. The energy interval e ~ Vq marks a transition from a "soft", small 
angle scattering to "hard" scattering, when propagation direction changes 
drastically in a single scattering event (in this case there is no need to dis- 
tinguish between / and /*). When the particle energy becomes smaller than 
Vq, and it is lowered towards the classical percolation threshold, reflection 
from potential barriers starts to play an important role and serves as "bot- 
tleneck" for transport. The transport slows down, because the particle can 
get temporarily trapped in a small region and it has to make many bounces 
between the potential barriers before it escapes. Eventually, at some critical 
energy Ec transport (at zero temperature) comes to complete halt. Clas- 
sically, such transition to an insulating state would have occurred at the 
percolation threshold Eper, with a diffusion coefficient [51] behaving in the 
threshold vicinity as D ~ Do[{e/ Eper) — 1]*, where t can be identified with 
the percolation conductivity exponent, close to 1.7 in 3(i [STl |52] . Classical 
treatment, however, is never adequate close to the transition, where quantum 
interference effects dominate. Thus, when the energy is lowered towards Ep^r, 
the diffusion coefficient, at first, will behave in accordance with the classical 
percolation theory but, eventually, a crossover to the quantum (Anderson) 
transition must take place 1^53] . This transition occurs at the mobility edge 
Ec, which is strictly larger than -Eper- (We stress that this consideration 
applies only if Vq >> Eq, whereas in the opposite case, considered in the 
previous subsection, no trace of classical percolation appears in the transport 
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picture.) The difference (E^ — Ep^r) can be roughly estimated by locating 
the energy at which the classical expression Do[{e/ Eper) — 1]* gets equal to 
the "diffusion quantum" h/m. Such estimate yields Ec = Eper{l + r/"^/^*). 

One can also appreciate the importance of quantum effects, while moving 
towards Ep^r from the low energy side. Well below Ep^r states are localized, 
essentially in a single potential well. When e is raised, a state spreads over 
several wells. At e = -Eper, a classical particle would get delocalized but 
quantum interference prevent delocalization. Only when e > Ec > Eper, does 
delocalization become possible. 

The above relation between Ec and Ep^r is based on the assumption of 
two well separated scales: a "microscopic" scale, at which local diffusion 
(with a diffusion constant Do) takes place, and a macroscopic one, related 
to the topology of the percolating cluster and responsible for the "slowing- 
down factor" [{e/Eper) — if- This assumption is by no means obvious even 
for the "generic" random potentials, usually assumed in the percolation the- 
ory |HT | 152 ) . Moreover, one can construct "unusual" potentials, the extreme 
example being a potential which, although random, is strictly zero within 
a connected spatial region. Percolation is then possible at all energies, and 
there is no connection at all between classical percolation and quantum lo- 
calization. Actually, a blue detuned speckle might be not so far from this 
extreme example since, due to the high probability for low values of the 
potential, the percolation threshold (in 3d) is only about 4 ■ 10~^Vo [54] . 

3.3 d = 2 and r/ « 1 

In two dimensions all states are localized. If, however, the disorder (at the 
given energy) is weak, i.e. kl » 1, the localization length is exponentially 
large. In this case, in a finite size sample, the wave function can easily 
spread over the entire sample, and the states can be considered as extended, 
for all practical purposes. When disordered increases (or the particle energy 
decreases, for fixed disorder), the localization length becomes smaller and, 
eventually, the regime of strong localization is reached. Thus, in 2d, one 
speaks about a crossover from weak to strong localization, instead of a strict 
transition between extended and localized states that is taking place in 3d. 

This picture can be rephrased in terms of transport. At high energies 
the ordinary diffusive transport is possible. Indeed, it would take unrealisti- 
cally large samples and exponentially long times (in addition to the almost 
complete absence of any inelastic processes) to find out that the initial wave 
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packet gets, in fact, localized. Under appropriate conditions, one can ob- 
serve weak localization corrections, on top of diffusion |55[ |33l |35]. When the 
energy is lowered, the diffusion coefficient decreases and one crosses over to 
the strong localization regime: the localization length becomes much smaller 
than the sample size and transport is inhibited. 

As in Sec. 3.1, we start with an estimation of the cross section for a single 
scatterer- a potential barrier of hight Vq and radius Rq. For a slow particle, 
kRo << 1, this cross section in the Born approximation is |42] a ~ r]'^/k and, 
unlike the 3d case, it is energy dependent. The approximation is valid for 
fci?o » rj^ ■ Next, we estimate the mean free path for a particle moving in 
the presence of many scatterers, in concentration riscat ~ -Rq ^- 

/ = :^ ^ (30) 

cr Tj^ 

Furthermore, the loffe-Regel criterion, kl » 1, results in the condition 
/ci?o » Tj (this is more restrictive than the condition kR^ » rf for a single 
scatterer). In terms of the particle energy one obtains [271 [56]: 

e » VqT] = E^. (31) 



This energy scale has already appeared in (26), although its significance here 
is quite different from that in Sec. 3.2, devoted to the smooth "id potential. 

Consider now the fast particles, kR^ » 1. The transport scattering 
cross section for those differs from a ~ r/^/A; by a product of two small 
factors (compare to the corresponding estimate in Sec. 3.1). The factor 
1/kRo accounts for the narrow scattering sector, in the forward direction, 
and the factor l/(/ci?o)^ accounts for the difference between the scattering 
and the transport cross sections. Thus, cr* ~ rf/k^R^, and the transport 



mean free path /* = Rl/a* ~ i?o(e/Vo)^, as in (17). 



The quantitative treatment of /*, in the Born approximation, is based on 



the 2d counterpart of (18): 



1 _ 1 

r* h 



Fo^z/^ (i(/)f |^2fcsin0 (1-COS0), (32) 



where v = m/2Txh^ is the 2d density of states and F is the Fourier transform 
of the correlation function T{R/Rq) for the 2d random potential. For a 
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Gaussian correlation function, F = exp(— i?^/i?Q), the result for /* = hkr* /m 
is (SUES EH : 

Here the function f{x) = e^[/o(x) — where Jq and Ji are the modified 
Bessel functions. For small x, f{x) ~ 1; for large x, f{x) ~ y/7r(2xY^'^. In 
these limits of slow and fast particles, respectively, one can immediately write 
down the corresponding values for the mean free path, in complete agreement 
with the above estimates. In terms of the diffusion coefficient, Z)e = vl* /2, 
the limiting values are: 



where, as in (20), we have made it explicit that the energy should be counted 
from the average potential value. 

The important case of a correlation function, corresponding to a speckle 
potential created by a uniformly illuminated circular diffusive plate, was con- 
sidered in [27J. Although that correlation function decays only as (with 
oscillations), the integral j T{Y{)(PR is finite, so that this case is, qualita- 



tively, not different from the Gaussian correlation. The equation (34) still 
holds, albeit with different numerical coefficients (note that our definition of 
Eq and 7] differs by a factor of 2 from that in [27]). 

Unlike the "id case, when the condition [e — V) » W was sufficient for 
diffusive transport, in 2d the condition (e — V^) >> E/\ is only a necessary 
one. In an ideal world of infinite samples, infinite times and zero tempera- 
ture, quantum interference effects would eventually take over and transport 
would cease. In reality, factors like dephasing or finite sample size suppress 
interference and provide room for diffusion. Interference effects show up 
in a correction to the diffusion coefficient (weak localization). As the just 
mentioned "ideal conditions" are approached, interference becomes more im- 
portant and the diffusion coefficient approaches zero. All this is discussed in 
great detail in [2^ ITT] . For {e — V) ^ a crossover to strong localization 
takes place. The localization length ^ rapidly decreases with energy, and for 
{V — e) » Ea the region of strongly localized states, well separated from 
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each other, is reached. This region of the spectrum has been aheady discussed 
in Sec. 3.1, one should only set c? = 2 in the corresponding expressions. 

In condensed matter physics, transport in the region of strongly localized 
states can occur only at finite temperatures, and it is due to hopping |3H 
an electron can hop from one localized state to another by absorbing a 
phonon. In principle, one should be able to observe a similar phenomenon 
for localized cold atoms, by modelling phonons with the help of a random, 
time-dependent optical potential (this remark applies, of course, also to 3d). 



3.4 d = 2 and 77 >> 1 

A few brief remarks will suffice to summarize the situation for this case. Using 



considerations entirely similar to those in Sec. 3.2, one can see that Eq. (28) 
for the transport mean free path is valid also in 2d, as long as e >> Vq. The 
corresponding diffusion coefficient is D = vl* /2 ~ Do(£/Vo)^''^, i-G- the same 
as in 3d, up to a numerical factor. 

When e approaches the typical barrier height Vq, the diffusion coefficient 
reaches a value ~ Dq, and it keeps decreasing when the energy is lowered 
further. The decrease occurs for two reasons: first, it becomes more a more 
difficult for a (classical) particle to find a path around the potential barriers 
and, second, the quantum interference effects are being enhanced, leading to 
strong localization. When D reaches a value of order h/m, the classical pic- 
ture breaks down completely. Diffusion of a classical particle in a 2d speckle 
potential has been numerically studied in |28j . and it would be interesting to 
account for quantum effects, including the classical-quantum crossover. 



4 A superfluid-insulator transition in a disor- 
dered BEC 

The behavior of a BEC in a random environment is a vast subject, with a 
long history. The question of how disorder can affect, and possibly destroy, 
superfluidity and superconductivity has been repeatedly discussed in the lit- 
erature [59| l33] . In particular, Ref. [59j contains a thorough discussion on 
^He absorbed in a porous Vycor glass. The porous, random structure of 
the medium can lead to a signiflcant reduction of the superfluid transition 
temperature Tg. 
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Ultracold atomic gases, which allow for a very accurate control over both 
the disorder and the interparticle interaction, provide an ideal system for 
studying the interplay between the two factors. The problem is of funda- 
mental interest, because it brings together two central themes in condensed 
matter physics: interactions and disorder. In the absence of interactions all 
bosons would drop into the lowest potential well, so the ideal Bose gas is 
clearly not the appropriate starting point for the theory. Any finite repulsive 
interaction has a drastic effect, pushing the particles apart and not allowing 
for too high local density n(r). If the average particle density, n, is finite, 
the particles will get spread over the entire (infinite) system. 

We will assume the model of a continuous random potential V(r) (see 
previous section), thus, leaving out the work on disordered (or quasiperi- 
odic) optical lattices and the disordered Bose-Hubbard model (see [5] for an 
extensive discussion and references, and [60kj61j for recent experiments). Fur- 
thermore, we will not discuss finite temperature effects or excitations above 
ground state [Ml ESI ESI El] but focus on the zero-temperature, quantum 
phase transition, i.e., the transition between a superfluid and an insulator 
under a change of the chemical potential /i (or the particle density n) |39lH0] . 
Ref. [40j contains a detailed exposition of various regimes, in 2d and 3c?, in- 
cluding the effect of a harmonic trap. We refer the reader to [40] for details, 
and here only review the transition in 3d, for a short range random potential, 
7] « 1. 

We consider, thus, N bosons in a box of volume Q, in the thermodynamic 
limit {N,Q — )■ oo,n = N/Q = const). Repulsive interactions and the "gas 
condition", na^ << 1, are assumed. Here a is the scattering length, related 
to the interaction constant as g = Airh'^a/m. In the absence of disorder, 
the chemical potential is fi = gn [T^ [T3] . A weak random potential can be 
treated in perturbation theory [65^, '66]. The important energy, characterizing 
a short range random potential, is W (see Eq.(16)). A particle with the de- 
Broglie wave length Ap >> Rq experiences not the "bare" potential V{r) but 
a smoothed, effective potential, and the scattering can be considered weak 
as long as Ap << / ~ Rq/t]'^. For a BEC the role of Ap is played by the 
healing length C,h = h/y/2jrign = h/ yJ2m^ and the effect of the potential is 
small, provided the condition « i-e. fi » W, is satisfied. Thus, for 
fi » W one has a weakly disordered BEC. It is described by a macroscopic 
wave function \E'(r) which, although fluctuating in space, maintains coherence 
over the entire box. The system is a superfluid. 

The perturbation theory is expected to break down when fi reaches a value 
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of order W, i.e. when particle concentration becomes of order W/g ~ 1/C'^a, 
where we have introduced the Larkin length defined in Sec. 3.1. At such 
concentrations disorder becomes more important, fluctuations grow stronger 
and at some point the coherence gets disrupted: the superfluidity is destroyed 
and the system undergoes a transition to the insulating state. Let's denote 
this critical concentration by ric and the corresponding value of the chemical 
potential /i = fic- Generally, there is no simple relation between fic and the 
mobility edge Ec (or the energy W). To determine fic one has to solve the 
many body problem of interacting disordered bosons, whereas Ec is a single 
particle property. In other words, there is no straightforward connection be- 
tween superfluidity and the nature of the single-particle states (localized vs 
extended). Indeed, it has been argued long ago [67] that, due to the interac- 
tion induced screening, a condensate can exist when the chemical potential 
is well below Ec- Nevertheless, an estimate of fic for weakly interacting dis- 
ordered bosons |39| HQ] shows that W is an important energy scale in the 
problem. The arguments and estimates in [3^ HO] were made for a Gaussian 
random potential. Since, however, the main "action" is taking place in a 
narrow interval near V, the considerations should apply also to a speckle 
potential. For small n, when /x (measured from the average potential value 
V) is negative and large in the absolute value, the distribution of particles 
in space is highly inhomogeneous. The particles, in the ground state, fill 
the potential wells with radius smaller than some value R (the correspond- 
ing depth is ~ —h^/mR^), and the purpose is to determine this value, for 
a given n. Let's denote by nw{R) the concentration of such wells. For R 
small, namely R << C, the typical distance between neighboring wells is 
exponentially large and their concentration is [40] 



where the preexponent / plays no significant role in the forthcoming estimate. 
The average number of particles per well is N^{R) = n/n^{R) and the 
particle density in a well is np{R) = 3iV^(i?)/47ri?^. The total energy of 
Nyj{R) bosons in such a well is comprised of two parts: the binding energy 
and the interaction energy. Introducing the energy per particle, E{R), one 
can write (with somewhat arbitrary numerical coefficients): 




(35) 



E{R) 



+ gnp{R) 



3h'^na 



(36) 



2mi?2 n^{R)mR^ ' 
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Now, using (35) and minimizing (36) with respect to R (for fixed n), one 
obtains [lO] : 

For this resuh to be vahd, one must require n « Uf.. On the other hand, 
since the number Ny^{R) must be large, n should be not too small, namely: 



n, 



:exp [-(£/3a)^/^] « n « ric, (38) 



which is consistent only if £ >> a. Under the condition (38), tunneling 
between the wells is negligible and the number of particles in each well is 
fixed. The state is insulating (a fragmented condensate). The chemical 
potential is estimated as 



/^H — ---Wln^p . (39) 

When n increases, the size R{n) of the relevant wells gets larger and, for 
n ~ rici it becomes comparable to the separation between the wells (both are 
of order C). Tunneling becomes significant and at some point a transition 



to the coherent superfluid state occurs. Of course, ric, as defined in (37), 
is not the precise transition point but only an order of magnitude estimate. 
Similarly, the critical value /Xc of the chemical potential is not accurately 
known. 



5 Dynamics of cold atoms in the presence of 
disorder 

The present section is devoted to the time evolution of atomic clouds sub- 
jected to a random potential. In a typical experimental set-up a ultracold 
atomic gas, of either bosons or fermions, is released from a harmonic trap, 
Vtr{T) = ma;^r^/2, and undergoes expansion while being scattered by a ran- 
dom potential V{r). The trap determines the initial state of the gas, whereas 
the potential ^(r) affects its subsequent time evolution. At some moment t 
an image of the atomic cloud is taken, providing information about the mode 
of transport. We discuss separately two atomic systems: BEC and a cold 
Fermi gas. 



26 



5.1 Ballistic expansion of a BEC 

We start with the problem of a free, balhstic expansion, i.e. the random 
potential is absent. BEC in the trap is described by a macroscopic wave 
function, satisfying the stationary Gross-Pitaevskii equation 

1 

- — + -mwV^^ + = /i^. (40) 

Usually, interaction among the atoms in the trap is of great importance: the 
size of the condensate, ~ {fi/mco'^Y^'^, is much larger than the "oscillator 
size" ao = {h/mcoY^'^. The latter is, roughly, the size of a single particle 
ground state wave function, and in the absence of interactions all the bosons 
would have condensed into this single particle state. The inequality i?^ >> flo 
implies » hu. This is the condition for validity of the Thomas- Fermi 
approximation which amounts to the neglect of the kinetic energy term in 



(40), thus, resulting in a simple expression (an inverted parabola) for the 



condensate density 

„(M=o) = i*(M=or = ^(i--;J) . R.={^y^ (41) 

where we have added the time argument t = 0, to emphasize that this is the 
condensate profile in the trap, prior to its release. The corresponding wave 
function is a square root of this expression, with an arbitrary overall phase 
which can be set to zero: 



^(r,t = 0) 



/ 2 \ 

/ r \ 



9 V rV\ 



1/2 



F(r). (42) 



This expression is accurate in the bulk of the condensate, and it breaks down 
only close to the boundary, for r ^ R^. 

Assume now that at t = the trap potential is switched off, and for 
t > the BEC undergoes free evolution, according to the time dependent 
Gross-Pitaevskii equation: 

ih— = A^ + ^I^P^ , dr\^\'^ = N. (43) 

ot 2m / 



with the initial condition given in (42). The wave function is normalized 



to the total number of particles A^. ^I/(r,t) is a complex function, \1/ 
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^/nexp{i6), i.e. the condensate acquires a dynamic phase 6{r,t). Eq. (43) 
can be rewritten as a set of two hydrodynamic equations, in terms of the 
condensate density n{r,t) and its velocity v(r,t) = (/i/m)V^(r, t) |T2l [T3] : 

dfi 

— + divn-v = (44) 
ot 

m— + V I -mf + gn j =0, (45) 

where in the second equation the "quantum pressure" term has been ne- 
glected, which is justified if the healing length is smaller than the charac- 
teristic scale over which the density is changing (see the discussion in [13] ). 
These equations have to be solved with the initial conditions v(r, t = 0) = 
and n(r, t = 0) given in (|4l|. The solution is [HI EHl i2| : 

"(r.<) = ^(l-j^). v(M) = ^r. (46) 

where the scaling factor b{t) obeys the ordinary differential equation b = 
u'^/b'^"^, with 6(0) = 1,6(0) = (a dot on b denotes time derivative). At 
= 2 its solution is b{t) = \/l + uH'^^ and at = 3 it is qualitatively similar, 
namely: 6 ~ 1 for t « 1/u and b ^ ut for t » 1/u (linear evolution). 

Thus, in the course of the expansion the condensate density ?T,(r, t) retains 
its shape of an inverted parabola whose size increases according to the scaling 
factor b{t). The important time scale in the expansion process is to = l/oj. 



For t « Iq the time evolution is dominated by the nonlinear term in (43) 
and the condensate energy is due primarily to the interaction. By the time 
if: ~ to most of the interaction energy gets converted into the kinetic energy 
of the condensate flow. The wave function \E'(r,to) exhibits rapid spatial 
oscillations which account for the large kinetic energy of the condensate. 



For t » to the interaction term (i.e. the nonlinear term in (43)) becomes 
negligible and the wave function evolves according to the linear Schrodinger 
equation. The function at t = to is, roughly, 

^(r, to) ~ fC-) exp(irVao) = ^{r), (47) 

where the argument r/2 in F indicates that by the time to the condensate's 
size has increased by a factor of 2 or so. One can check that the kinetic 
energy in this wave packet, E^in = |~ / |V\l'p(ir, is of the same order of 
magnitude as the interaction energy, Eint = f / l^l^c?!"- 
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5.2 Diffusion of a BEC 



Here we study the time evolution of a BEC, in a random potential V{r), upon 
its release from the trap. The evolution is governed by the Gross-Pitaevskii 
equation, with the random potential: 

ih— = /\^! + V{r)m + g\m\^m (48) 

Sometimes it is possible to disentangle the effects of nonlinearity and dis- 
order [TOl [5? j. For instance, one can switch the trap off at t = 0, let the 
condensate freely expand for a time interval of the order of to = l/i^? and 
only then switch on the random potential. In this way the evolution is sep- 
arated into two stages. The first stage, for times t ~ to? is dominated by the 
nonlinearity (there is no disorder yet), whereas at the second stage, t » to, 
the nonlinearity can be neglected and the evolution proceeds according to 
the Schrodinger equation 

= + ^W"^ ' > ^o)' (49) 



with the initial condition \I'(r,t = to) = *&(^), see (47). Setting the start of 
the second stage at t = to is, of course, only a rough estimate (one should 
wait for a time few times larger than to before switching on V^(r)) but it does 
capture the main idea of the two-stage evolution. 

In a different set-up, one releases the condensate, at t = 0, directly into a 
random potential (rather than switching on the potential later) and the time 



evolution for t > follows the equation (48). If V{r) is sufficiently weak, 
one can still argue in favor of a two-stage scenario. During the time t ~ to 
nonlinearity dominates over disorder and the expansion is essentially ballistic 
(Sec. 5.1). For t >> to the nonlinearity is weak and disorder becomes the 



important factor: the evolution is close to linear, according to (49). Whether 
a given potential V^(r) is weak or strong depends on the wave number k of 
the relevant component of the condensate wave function. If, for instance, the 
chemical potential in the trap is much larger than the mobility edge Ec, 
then most of the atoms, upon their release from the trap, will experience weak 
disorder, because the parameter k^l is large (A;^ = \j2m^/h'^). On the other 
hand, spreading of the small-A; components might be inhibited by disorder 
and they will stay localized in the vicinity of the trap. For these components 
the nonlinearity might remain important for all times. In any case, the first 
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set-up, i.e., when the potential is switched on after sufficiently long interval 
of a ballistic expansion, is better suited for the two-stage scenario. 

Assuming the validity of this scenario, it remains to study the linear evo- 
lution (49 ), with the initial condition (47). This condition carries information 
about the initial state of the condensate in the trap and about the first stage 
of expansion, dominated by the nonlinearity. To simplify the notations, we 
reset the time to to zero, so that the formal solution of (49) is 

^(r,t) = y"rfRG'(r,R,t)$(R), (50) 

where G is the retarded Green's function of the Schrodinger equation (49). 
The average particle density 

n(r,t) = |^(r,t)|2 = JdnJ dR'G*{r, R, t)G{r, R', t)$*(R)$(R'). (51) 

The product of the Green's functions can be Fourier transformed as 



G*(r,R,t)G(r,R',t) = J^J ^e-^G*(r, R, e + ir])G(r, R', ^ - i^]) . 

(52) 

The product in (52), for weak disorder, can be calculated in the diagram 
technique [35j : 



G*(r, R,e + ^Q)G{r, R', e - ^Q) = -2P,(r, R, n)ImG{R - R' , e), (53) 
where the diffusion ladder is defined as 



P,(r, R, n) = -^G*{r, R, e + ^fi)G(r, R, e - ^Q) (54) 

and the average Green's function 

G{p,e)=Go{p,e)e-^, (p = R-R'). (55) 



Here Go is the free Green's function and and are, respectively, the 
(average) density of states and the mean free path. Putting all the pieces 
together, we arrive at the following expression for the particle density: 



n{r,t) = -^ f dR f dp f rfeP,(r, R, t)/mG(p, e)<l>*(R + |)$(R - |^ 



(56) 
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where the quantum diffusion kernel Pe(r,R, t) is the Fourier transform of 
the diffusion ladder P(;(r,R, fi). Although this result was derived for weak 
disorder, scaling arguments (as presented e.g. in |33t I7T]) suggest that, with 
the appropriate form for the kernel P, the result is valid also when disorder 
is strong. Transforming G{p,e) to the momentum representation, G(k, e), 
enables us to write 

n(r,t) = J dnj dkj deP,{r,R,t)ImG{k,e)W{k,R), (57) 
where the Wigner function is defined as 

W{k,R) = ^,j c?pe^''-''$*(R+ f )$(R- f ). (58) 



Eq. (57) contains three ingredients: (i)The Wigner function carries infor- 
mation about the initial condition. (ii)The factor — (l/7r)/m(j'(k, e) = v4(k, e) 
is the spectral function. For a free particle it is equal to 5{e — ek), with 
ek = h'^k'^/2m (the "on-shell" relation between the energy and momentum). 
Disorder broadens the (5-function into a Lorentzian or, if sufficiently strong, 
into some more complicated shape, (iii) The quantum diffusion kernel prop- 
agates a particle, with energy e, from R to r, in time t. 



Eq. (57) simplifies in certain limits. If one uses 5{e — ek) for the spectral 



function, which can be justified for sufficiently weak disorder, one obtains 
n(r,t) = J dnj rfkPfc(r,R,t)W^(k,R), (59) 



where = P^^. . This equation appears, in a somewhat different form, in [2] 
and it has a simple interpretation: The initial distribution in phase space, 
Ty(k, R), is propagated in time by the propagation kernel Pfc(r,R, t). One 
can rewrite this equation in terms of the momentum p = hk, making it look 
entirely classical. Note, however, that both W and P^ are quantum objects. 

Another type of approximation, appropriate for large distances r, is based 
on neglecting the P-dependence of Pfc(r, R, t). For particles which had prop- 
agated a distance r » R^, the trap (i.e. the initial size of the cloud P^) can 
be viewed as a "point object" and P can be set equal to zero (the center of 
the trap). The integral of W{k, R) over R gives the momentum distribution 
|0(k)p, where 0(k) is the Fourier transform of the initial wave function $(r). 
(57) now reduces to 

n{r,t) = J j^jm' J deP,{r,0,t)A{k,e). (60) 
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A still simpler expression for n(r, t) is obtained by combining the two afore- 



mentioned approximation, i.e. by equating the spectral function in (60) with 
5(e-ek): 

— \m'Pkir,t), (61) 
where the second argument, R = 0, in Pfc(r, R, t) is omitted. The expression 



(61 ) can be visualized as dynamical evolution of a cloud of classical particles, 
which are initially at the origin and whose velocity distribution is determined 
by |(/)(k)p. Thus, (61 ) applies either to a BEC (a single coherent wave packet) 



or to a cloud of independent, non-interacting particles, e.g., a Fermi gas (see 
Sec. 5.4). The information about the system is in the initial distribution 
|0(k)p, which can be, for instance, the thermal distribution (Bose or Fermi) 
at a given temperature. 



The interpretation of |\E'(r, t)]"^ in terms of "particles" is not limited to (61 ) 
but holds generally. One can imagine performing integration over R and k 
in ( [57| ), ending up with the spectral decomposition of the condensate density 
n{r,t) = f |^ne(r,t). Similarly, integration over R and e would produce the 
momentum decomposition J ^^^?2k(r,t). Thus, one can interpret the total 
density n(r, t) as being composed of fictitious particles, or atoms, with a 
given energy (or momentum), and in what follows we often use this intuitive 
picture, talking for instance about fast or slow particles. One should keep in 
mind, though, that the condensate wave function is a coherent unit and that 
only the total density has the straightforward interpretation in terms of the 
actual bosonic atoms. 

To proceed further one needs an expression for the quantum diffusion 
kernel Pe(r, R, t). This is the main object of the theory because it determines 
the dynamics, for a given e. For weak disorder, when propagation is by simple 
diffusion, P is the standard diffusion propagator: 

where is the (energy dependent) diffusion coefficient (see Sec. 3). When 



e approaches the mobility edge Ec, the disorder gets stronger and (62) is 
not applicable any longer. For that case one can use the scaling theory for 
the Fourier transformed kernel, Pe(r,R, fi), which is related to the dynamic 
diffusion coefficient D^{fl) [33j . 
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We start with the ordinary diffusion, (62), and restrict ourselves to the 



long time limit, when (61) is valid. If D were constant, i.e. energy indepen- 



dent, (61) would reduce to a Gaussian function (the momentum distribution 



is normalized to the total number of particles J 



dk 

(27r)" 



N). Due to 



the energy dependence of D, the shape n(r, t) deviates from Gaussian and 
it depends on the initial momentum distribution |0(k)p. The characteristic 
width of this distribution is fc^ = ^j7im\ifW = and it is due to the 
rapidly oscillating part which the wave function acquires after the first stage 
of a ballistic expansion. To a good approximation, in any space dimension, 
the fc-dependence of |0(k)p is determined by the factor (1 — ^), 
For instance, in 2d 



Q{k^V2-k). 



(63) 



Let us consider a short-range correlated potential, i.e. rj << 1, and 
assume that fc^-Ro << 1- It must be also assumed that the chemical potential 
/i of the condensate (prior to the release from the trap) is much above the 
mobility edge Ec, so that the great majority of the fc-components experience 
weak disorder and propagate by diffusion. In 2d the stated conditions are the 



same as given in the inequalities in the first line of (34). The mean free time 
is then independent of k and the diffusion coefficient is = h^k'^r /2m?' = 
D^{k/k^y', where D,, is the diffusion coefficient for k = k,, 



The upper limit in the integral (61) is k^\/2. Since for sufficiently small 
k diffusion breaks down, a lower cutoff, kc, should be introduced in the 
integral. The cutoff is estimated from the condition kl = 1 which, using 
(30), gives kc = {ti/Rq) = {m/fiTY^'^ . This cutoff is a rough benchmark 



between the components which diffuse away and those which stay localized. 

tegral i 
^c, yields 



Writing the kernel (62 



with the expression (63) for |0(k 



in terms of k and calculating the integral in (61) 
P and with a lower limit fc, 



n r 



N 



[l + a) 



EAa)-Ei[ 



,2a. 



+ 



-5^ exp( 



2a 



(64) 

where Ei{x) is the exponential integral, a = r'^/8D^t and 6 = {k^/k^) = 
^yh/2jj,T is the dimensionless cutoff. The cutoff is important only for small 
r, to avoid the divergence at r = 0. For >> ht/m the parameter 2a/5^ is 



large and the two 5-dependent terms in (64) can be neglected. 
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In Fig.2 we plot the normalized particle density, n{r] 
for some fixed time t, as a function of the normalized distance f 



47rD^tn{r,t)/N, 
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Figure 2: Density as a function of position (both quantities are normalized 
as explained in the text). The cutoff 5 is zero for the solid curve and it is 
0.04 for the dashed curve. 



The dashed curve corresponds to a cutoff 5 = 0.04 while for the solid curve 
6 = 0. The curves exhibit an exponential tail (with a power-law prefactor) for 
large f. This tail comes from the large-A; part of |0(k)p, i.e. from particles 
which, in a given time, manage to diffuse far away. It can be obtained 
analytically, using the large-^ asymptotics Ei{z) ~ exp{—z). The sharp 
increase in h{r) for small f is due to the small-Zc components whose diffusion 
coefficient decreases as k"^. Since for small z, Ei{z) ~ — Inz, the solid curve 
has a logarithmic divergence at the origin. The divergence is rounded off by 
the finite cutoff 6. 

It is sometimes convenient to replace the sharp cutoff in ( [63| ) by a smoother, 
Gaussian cutoff, thus, taking |0(k)p = (47rA^/A;^) exp(— A;^/A!^), with the re- 
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suit [ST]: 




where Kn is the zeroth order modified Hankel function, and no lower cutoff 



5 was introduced. This result is very similar to (63), with 6 = 0. Indeed, the 
function Kq{z) has an exponential tail (with a power-law prefactor) for large 
z, and it diverges logarithmically for small z. 

It is interesting that is the only system-dependent quantity (except 



for the total number of atoms N) which appears in (64) or (65). Using 
the natural quantum unit of diffusion, h/m, one can define a dimensionless 
parameter mD^/h ~ HT*/h, where r* is the transport mean free time for a 
particle with energy fi (this relation is not restricted to the above 2d example). 
This parameter carries information about the initial state of the condensate 
(/i depends on the nonlinearity and on other factors), as well as about the 
subsequent dynamics (via r*). It is analogous to the parameter Epr/h which 
determines the transport properties of a disordered metal {Ep is the Fermi 
energy) . 

Diffusion of cold Kb atoms, in a 2(i-geometry, was observed experimen- 
tally in [9j. Motion in the vertical {z) direction was confined and transport 
in the x, y-plain was studied. The temperature was somewhat above the 



condensation temperature, so that instead of the inverted parabola ( 63 ) one 
should use the appropriate thermal distribution [9j. Since the shape of the 
diffusing cloud is not very sensitive to the precise initial distribution, pro- 
vided it still contains fast and slow particles, one can expect results similar 
to those in Fig.2. Indeed, profiles of the (integrated along x or y) density, 
shown in Fig. 3, have clear resemblance with Fig.2. Note the difference be- 
tween curves in (b) and (c) of Fig. 3, for the same propagation time. It is due 
to the anisotropy in the potential and hence in the diffusion coefficient. The 
experimental data in [9] were fit numerically to a model with a anisotropic, 
energy-dependent diffusion coefficient. 



35 
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Figure 3: Atomic column density after planar expansion of an ultracold 
gas in an anisotropic speckle potential, (a) Image after 50 ms of expan- 
sion, (b), (c) Integrated density along the two major axes. The plain dots 
(open squares) correspond to 50 ms (200 ms) of expansion. (Reprinted, with 
permission, from Phys. Rev. Lett. 104, 220602 (2010)). 

The characterization of a diffusing atomic cloud is not limited to the aver- 
age particle density n{r,t). Indeed, in a given realization of randomness the 
density pattern n(r, t) exhibits spatial fluctuations, which are smoothed out 
only in the process of averaging over many realizations. The fluctuations are 
due to multiple scattering of the matter wave on the random potential. The 
phenomenon is similar to optical speckles for the electromagnetic waves (see 
Sec. 2) and can be termed "matter wave speckles". The density fluctuations 
can be described by the (equal time) two point density correlation function, 
n{r,t)n{r',t), as well as by higher order correlation functions. The function 
n{r,t)n{r',t) reveals both short and long-range correlations, which ( for a 
white noise potential) have been studied, respectively, in |72) and |73) . 

Finally, let us recall that our discussion has been limited to the "two-stage 
scenario": first, a free ballistic expansion and only later, when the nonlin- 
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earity had decreased sufficiently , the disorder is switched on. An interesting 
and conceptually important question is to what extent the diffusion is mod- 
ified in the simultaneous presence of disorder and nonlinearity, i.e. when a 
wave packet evolves according to (48). This problem has been addressed in 
[7H [THl [70]. It was shown [73] that in the first order in perturbation theory 
the nonlinearity causes a slight renormalization of the diffusion coefficient 
which controls the condensate spreading. In the works [75l 176] , which went 
far beyond the simple perturbation theory, a generalized diffusion equation 
for the spectrally resolved (average) density n^{r,t) was derived. The "dif- 
fusion coefficient" in that equation depends on r and t, through the total 
particle density n{r,t), which has to be determined self-consistently. 



5.3 Localization of a BEC 



The general expressions (56), (57), derived in the previous subsection, apply 



also for strong disorder (i.e. when e is in the vicinity, or below the mobility 
edge Ec), provided one uses the appropriate expressions for the propagation 
kernel and the spectral function. It is convenient to define the Fourier trans- 
form, Pe(r, Q), of the kernel Pe{r, 0, t) = Pe{r, t). For ordinary diffusion, (62 ), 
Pe(r,f2) = (AnDf^r)'^ exp(— r^— ifi/D^). More generally. 



A7TD,{n)r 



exp 



-in 



(66) 



where D^{Q) is the dynamic, frequency dependent diffusion coefficient which 
carries information about transport beyond ordinary diffusion. This quantity 
has been extensively studied in the context of disordered electronic systems, 
where it is proportional to the ac conductivity 0"e(f2), with e being the Fermi 
energy [33 ] 171 ) . The i7 — > limit corresponds to the dc value. In the localized 
regime, e < E^, this value is zero, D^{0) = = 0. The dynamic diffusion 
coefficient, however, is not zero and, for small Q, it is given by 



{e < E, 



(67) 



where ^ is the localization length. Close to the transition ,^(e) oc {Ec — e)~^, 
where the exponent u ~ 1.57 j33]. The expression ( [67j ) is easily understood in 
terms of a^{Q): although there is no dc transport (the eigenstates at energy e 
are localized), an external field at frequency Q polarizes the medium, causing 
an oscillating polarization current. 
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Note that (67) is valid only for Q < ^dO^ where the crossover frequency 
flciO^ estimated below, decreases to zero when ^ — )■ oo. At the transition 
point, e = Ec, (67) has no region of applicability and it is replaced by Dc = 
(—i'j^lY^^, where 7 is a model-dependent coefficient. This behavior follows 
from the scaling theory of localization [771 [33] and it implies that at 
the mobility edge a wave packet spreads according to oc t^^^ (anomalous 
diffusion), rather than oc t ( ordinary diffusion). Comparing Dc with D^fl) 
in (67) yields fic(0 ~ \/l/^^-> ^ characteristic time tc = (l/Qc) This 
is the crossover time between anomalous diffusion and localization: A particle 
(a wave packet) with e somewhat below will propagate by anomalous 
diffusion up to time of order tc and only then (i.e. after it had spread over a 
region of size ^) will it "realize" that it is in fact localized. 

A similar crossover occurs at the other side of the transition, e > Ec, 
where the eigenstates are extended. Close to the transition there exists a 
macroscopic length ,^(e) oc (e — Ec)~'^, with the meaning of being a crossover 
length between anomalous and ordinary diffusion. A particle with an energy 
slightly above Ec first experiences anomalous diffusion, up to a scale ~ ^, and 
only at a larger scale exhibits ordinary diffusion, with a small dc diffusion 
coefficient = 1/hvc^ oc (e — E^'^ (here Vc is the density of states at the 
energy Ef)). Thus, slightly above the transition one can view D^Vt) as a sum 
of two contributions: the dc part and the anomalous part, proportional to 
{—iVty/'^. On the other hand, well above Ec-, where the disorder is weak, 
the diffusion coefficient (up to small fi-dependent corrections, responsible 
for the weak localization effects) is given by the expressions (20) or (29), 
depending on the value of the parameter rj. The existence of a macroscopic 
length ^(e) oc (|e — Ec\^^ on both sides of the transition is a consequence of 
scaling and it is quite natural: only by observing the system at a scale larger 
than ^ can one distinguish between localized and extended states. Using the 
expressions for D^VL) in various limits, and possibly extrapolating between 
them, one can reconstruct, via (66), the quantum diffusion kernel P^r^t). A 
useful interpolation scheme is the self-consistent theory j3H], employed in the 
numerical part of jH]. Its drawback, though, is that it yields a wrong value 
for the exponent v (1 instead of approximately 1.57, in 3d). 

In order to use (60) one has to know the spectral function, in addition to 
the quantum diffusion kernel. The average Green's function, in the momen- 
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turn representation, can be written as [33] (the Dyson equation): 



where e) is the self-energy. The simplest approximation for the self- en- 

ergy is the first order Born approximation: YIO^^ ^) — I (27r)~'^(i^/c'r(k — 
k')Go(k',e), where r(q) is the Fourier transform of the potential correla- 
tion function r(R/i?o) and G'o(k, e) is the unperturbed (i.e. in the absence 
of disorder) Green's function. A more elaborate approximation is the self- 
consistent Born approximation, which amounts to replacing G'o(k, e) by the 



full Green's function, (68), thus, obtaining a closed equation for the self- 
energy: 

(2vr)3e-ek.-E(k',e)- ^^^^ 

The solution of this equation is, obviously, an even function of Vq. Therefore 
the self-consistent Born approximation is hardly appropriate for the speckle 
potential, where the odd powers of Vq are known to be generally important 



EC'. 



An analytic solution of (69) can be obtained in the limit of a 
short range potential, when V^^r(R/i?o) can be replaced with u5(R), so that 
r(k — k') = u. The constant u is proportional to VqU^ and the formal 
limit Vq — 7- 0,-Ro ^ oo? Kf-Rn ~ const is implied. For instance, for F = 



-'^o 

exp{-R^/Rl) (see Sec. 3.1), u = tt^/ttV^RI The solution of (|69j), with 
F(k — k') = u, yields a k-independent self-energy 

The value of e* depends on the ultraviolet cutoff, needed to regularize the 



integral in (69 ). For small but finite correlation radius Rq the natural cutoff is 
1/Ro, which, in the leading order in the small parameter t], yields e* ~ —f]Vo. 
The mean free path (in 3d) is / = nh'^/mu'^ and it does not depend on e. This 
expression for / is entirely consistent with the result / = ARQ/rfy/TT of Sec. 



3.1 (see the line below (19)), if one substitutes the value u = 7i^^/7^VQRl. 
The mean free time is related to / in the usual way, = ml/hk, with 
hk = ^Jlmie — e*). From the condition kl = \ one obtains for the mobility 
edge Ec = e* ^ {h^ /2mP). For r] small, is negative (see the corresponding 
remark in Sec. 3.1). In what follows we set e* = (the edge of the spectrum) 
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and count all the energies from this point. The mobility edge becomes 



Knowledge of the self-energy enables us to write down the expression (68) 
for G{]<., e), and hence for the spectral function 



^/mG(k,e)^A(k,e) = ^- , , , x = ^, (71) 



7ri?c(x-i-|^)2 + x E, 



c 



where, instead of e, the dimensionless variable x has been introduced. 

Let us fix some distant point r and observe the atoms arriving at this 



point in the course of time, according to (60). If the mobility edge E^. is 
well below the chemical potential fi (recall that information about /i in (60) 
is carried by the momentum distribution 0(k)), then n{r,t), as a function 
of time, will exhibit several regimes. First, the fastest atoms will arrive, by 
ordinary diffusion (Sec. 5.2) and the density will reach a maximum at a time 
tarrivai — ('"^/S-D^). Later, slower particles, with e closer to Ec, will show up. 
Still later, the anomalously diffusing particles will make their appearance and, 
finally, for times much larger than the characteristic time tioc — {hr^/EJ^) 
only the localized part of the condensate will remain. All this has been 
studied in some detail in jH]. We discuss only the t — )■ oo limit, when only 



the localized part of the condensate is left, at any finite r. In this limit Df:{Q) 



is given by (67) which, upon substitution into (66) and Fourier transforming 



with respect to f2, yields the kernel 



P,(r, R, t ^ oo) = ^^^^^r^ (-^^ ) • (72) 



The density n(r,t — t- oo) = nioc{r) is then obtained from (57) by cutting the 
upper limit of the integral over e at e = E^. Since close to Ec the localization 
length ,^(e) is large, the localized part of the condensate exhibits long tails. 
For r well away from the initial location of the condensate, and making use 



of (71), one obtains |44J: 



niUr)-f{Ec/f,)-.i-] , (73) 



3 \ 



where f{z) ~ z^^"^ for z « 1 and f{z) approaches a constant for z >> 1. 
This localized tail contains information about the localization length expo- 
nent u. 
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An important quantity is the (average) number of locahzed atoms, A^/oc, 
or the fraction of the condensate, Nioc/N, which remains locahzed after the 



mobile part had diffused away. The expression for Nioc is derived from (57), 
by integrating it over r and setting the upper limit in the integral over e at 
Ec- Since integration of the kernel Pf^{r,t) over r gives unity, one obtains 



dk 



Nioc = J j^jm'J deA{k,e). (74) 

In view of the importance of this equation we give an independent, and rather 
general, derivation. Denoting by e^, \a > the energy eigenvalues and eigen- 
vectors for a given realization of disorder, one can write \l/(r,t) =< r|\E' >= 
< r\a >< «!$ > exp{—ieat) , where $ is the initial wave function. In the 
t — )■ oo limit, for any fixed value of r, only localized eigenstates contribute to 
this expression. Taking its square modulus, integrating over r and using the 
orthonormality of the a-basis, yields a simple, intuitively obvious expression 
for the number of localized atoms Nioc = Yl' I < '^1'^ > where the prime 
indicates that the summation is over localized states only. This expression 
can be written as 

pEc _^ rEc 

Nioc = / dey^ 6{e-ea) < >< >= / de < ^\ImG{e)\^ >, 



TT 



— oo 



(75) 

where G{e) = (e — H) ^ is the resolvent operator (for the given disorder 



realization). Writing (75) in the momentum representation and averaging 
over the disorder results in Eq. (74). |0(k)p in this equation is a function 



of k/kfj_. It assumes a constant value ~ N/{k^Y ^r k « k^^ with a cutoff 
at A; ~ k^. Furthermore, A(k, e), as a function of A;, decays on a scale 



k^ ~ y/2mEjh (see (71)). It follows then from (74) that Nioc is a function 
of the parameter k^/kc (or ^/Ec). We will not compute the detailed shape 
of this function but rather analyze the two limits: 

(i) kfj_ » kc- In this case the effective region of integration is A; ~ kc, 
so that the function |0(k)p can be replaced by its value at k = 0, and the 
remaining integral over k gives the density of states z/(e), i.e. 

Nloc = m = 0)P dev{e) ~ N{kc/k,f. (76) 

^ — oo 

The fraction of localized atoms Nioc/N << 1. 
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(ii) fc^ << kc- Now the spectral function is "flat" in the effective integra- 
tion region, k < k^, and k in A(k, e) can be set to 0. The remaining integral 



over k in (74) is equal to A^, by normalization, and using (71), one obtains: 



loc 



N 



loc 



TX 



dx- 



X 



0.45. 



(77) 



This result is surprising since the natural expectation is that, in the limit 
of strong disorder, the fraction of localized atoms should approach unity. A 
short discussion of this result is given in The issue deserves further 



investigation, especially since ( 71 ) becomes doubtful in the strong disorder 
limit. 

Recently an experimental observation of Anderson localization of a BEC 
in "id has been reported jTT]. A dilute Rh condensate was released from 
an optical trap and allowed to freely expand for a time to = 50ms. After 
that time the random speckle potential was suddenly switched on and the 
dynamics of the condensate in that potential has been studied, i.e. images of 
the atomic cloud at various times have been taken. In Fig. 4 two sequences of 
images are shown (the time is measured from the instance when the disorder 
was switched on). The upper sequence corresponds to weak disorder, when 
a large fraction of the condensate evolves by diffusion. In the lower sequence 
disorder is strong and a significant fraction of atoms remains localized. 



l/R/h= 135 Hz 
S * 




Figure 4: 3(i expansion of the atomic cloud for the two values of the disorder 
amplitude. For small disorder (Vo//i = l?)bHz) one observes an essentially 
diffusive time evolution. For larger disorder (Vo//i = 680if2;) a significant 
fraction of atoms remains localized. (Reprinted from the arXiv: 1108.0137.) 
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The theory presented in [TT] starts with the following expression for the 
average density of atoms: 



n r 



t) = JdTiJ rfePe(r,R,t)D(R,e) , (78) 



where the function D(R, e) is interpreted as the semiclassical density in the 
position-energy space. Note that with D(R, e) = / (ikl4^(k, R)A(k, e) Eq. 



(78) is compatible with (57). For strong disorder, Ec >> /i, the spectral 
function is a broad, slowly decaying function of k, and it can be replaced 
by its fc-zero value, A{k = 0, e) = /(e). The remaining integral over k gives 
the initial (i.e. at time to ) density no(r), so that D(R, e) = no(r)/(e). A 
further approximation made in [H] was to replace Pe(r,R, i) by 6{r — R), 
thus, obtaining for the localized part of the condensate 

rEc 

fiiocir) = rao(r) / def{e) = no{r)fioc , (79) 

J —oo 

which is just the initial density profile, rescaled by the factor fioc- The func- 
tion f{E) was calculated in [llj numerically by diagonalizing the Hamilto- 
nian, for the realistic anisotropic speckle potential, and averaging over several 
realizations. A value fi^c ~ 0.22 was obtained, in the strong disorder limit 



(compare to 0.45 in (77), for the white noise disorder), and after some heuris- 
tic energy shift a fair agreement with the experimental results was achieved. 
Quantitative theoretical results are quite sensitive to various approximations 
made and, in particular, we note that replacing Pe(r, R, t) by 5(r — R), which 
might be reasonable for strongly localized states, will break down close to Ec 



where the localization length becomes large. Therefore (79) for the localized 



density profile cannot be the full story: n(r, t — )• oo) should have power law 



decaying tails (see (73)), associated with energies close to E^. 

Let us also mention a different type of localization, which takes place in 
the momentum space, rather than in real space. When a quantum system 
(chaotic in the classical limit) is "kicked" by a sequence of short pulses, the 
expectation value of the momentum square, < p^(t) >, evolves in time. In 
the well studied kicked rotor model < p^(t) > saturates in the long time limit 
[78j- a phenomenon known as dynamical localization, which was observed ex- 
perimentally in a system of cold atoms |79| . A 3d extension of the model 
[80] exhibits a richer behavior, analogous to the Anderson transition: De- 
pending on the kicking strength, K, < p^{t) > can either linearly grow with 
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time (diffusion), saturate in the long time limit (dynamical localization) or 
display the above mentioned anomalous diffusion, < p^(t) > oc t^^^ (critical 
behavior at the Anderson transition point). The latter occurs at a critical 
value, of the kicking strength. Such "dynamical Anderson transition" 
has been experimentally observed in [STl IS3]- A system of cold cesium 
atoms was acted upon by a pulsed laser wave. The authors have identified 
the transition point Kc-, measured the localization parameter, {Kc — K)~'^, 
close to the transition, and extracted the value of the exponent u, in very 
good agreement with the accepted numerical value. 



5.4 Cold Fermi gas 

Cold Fermi gases are also of interest, as reviewed in |Mt [55] . Although in- 
teraction between the fermions are generally important, in some cases they 
have only a minor effect on the dynamics of the atomic cloud. This is partic- 
ularly true for a polarized Fermi gas when the Pauli principle suppresses the 
main mechanism (the s-scattering) for the interaction. In this subsection we 
briefly review some phenomena pertaining to the expansion of a cold Fermi 
gas released from a harmonic trap. 

In the trap, the fermions occupy the single-particle states, 0„(R), of the 
harmonic oscillator potential The derivation of the expression 

for the average particle density n(r,t) proceeds as in Sec. 5.2. The only 
difference is that, instead of the condensate wave function $(R), we now 
have N incoherent single-particle contributions. At zero temperature, when 
the fermions occupy states up to the Fermi energy Ep, the corresponding 
Wigner function is 

Vr^(k, R) = -^0 (ef-^- \muj''RA , (80) 



{2'kY \ 2m 
and the dynamics of the gas, upon its release from the trap, is described by 



(57), with WV(k, R) instead of 14^(k, R). In particular, for weak disorder, 
when the spectral function can be approximated by 6{e — e^), one obtains 
|86j 



which is the precise analog of ( 59 ) . Diffusion of a Fermi cloud in a 2d speckle 



potential was analyzed in [86]. The authors used (81) with the diffusion 



44 



propagator (62) and with appropriate for the 2d speckle |27| . The den- 
sity profile f2(r, t) turns out to be very similar to that for a diffusing BEC 
(compare (64) of Sec. 5.2 with Eq.(26) in [86]). This coincidence is not sur- 
prising because the Wigner function ( 80 ) for fermions has much in common 



with that for a BEC . Indeed, integrating (|80j) over R produces, in 2(i, the 
momentum distribution |0(kp = (ao/47r)(fc|, — A;^), where k\ = 2mEp/fv^ 
and ao = {h/muY^'^ is the "oscillator size". Similarly, integration over k 
results in an inverted parabola {1/ 4:TTaQ){B?p — E?) for the initial density of 
the cloud {R\ = 2Ep/mu)'^). Thus, both the density and the momentum 
distributions have the same shape as for a BEC in a harmonic trap, see, re- 
spectively, (41) and (63). In 3d the shapes, although somewhat different, are 
qualitatively similar. The treatment in j86] can be easily extended to finite 
temperature T if one uses in (81 ), instead of the step function, the expression 
{exp[l3{^ + \muj'^R^ - + where (3 = l/fc^T. 

The problem of localization, as reviewed in Sec. 5.3, should apply, with 
slight modifications, also to the Fermi gas. So far, however, no specific cal- 
culations along these lines have been done. On the experimental side, 3(i 
localization of a Fermi gas has been observed in [TO], where a gas of ^^K was 
released from a harmonic trap into a speckle potential. According to the au- 
thors of [10] their speckle potential is well described by a anisotropic Gaussian 
autocorrelation function (rather than by the more complicated function in 
Eq. ([s])). The gas was spin polarized, in order to make the interatomic inter- 
actions negligible. The temperature range was between 200 and 1500nK, well 
above the Fermi energy, so that quantum degeneracy played no role. Upon 
release from the trap, some part of the atomic cloud had diffused away and 
the remaining localized component was observed. Again, as for the BEC |11] . 
the present state of the theory of Anderson localization in a anisotropic "id 
speckle potential does not allow for a detail comparison with the experiment. 

It is worthwhile to emphasize that the notion of the average particle den- 
sity n(r, t) involves two kinds of averaging. First, one performs the quantum 
(or thermal) average, i.e. writes down the expectation value (n(r,t)) of the 
density operator n(r, t), for a given realization of the random potential. Then 
averaging over the statistical ensemble of realizations is performed, yielding 
{n{r,t)) = n{r,t) [87] (In [86] this quantity is denoted as n{r,t), without 
the overbar). One should keep in mind, as has been particularly emphasized 
in |88| 189] , that in a single imaging experiment (with sufficient resolution) 
one does not measure the expectation value {h{r,t)) but rather one partic- 
ular event, i.e. some particular density pattern, n{r,t), whose probability 
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is dictated by the many-body wave function of the system. Therefore even 
in the absence of disorder a single experimental image will look "noisy" and 
grainy. (Only upon averaging over many measurements, taken under identi- 
cal experimental conditions, will one obtain the quantum expectation value 
(n(r,t)).) Fluctuations and correlations in such noisy density patterns are 
characterized by the (equal time) correlation function {h{r,t)n{r' ,t)) . For 
a clean, homogeneous Fermi gas, in equilibrium, this correlation function 
is described, for instance, in [90]. It exhibits decaying oscillations, with a 
characteristic spatial period Axq ~ kp^. Such oscillatory behavior occurs 
also in a gas confined to a harmonic trap. When the gas is released from 
the trap, it expands ballistically (in the absence of disorder). The size of 
the atomic cloud grows linearly with time and so does the correlation length 
Ax{t) ~ utAxQ. (In Id, an exact expression for the correlation function 
of a freely expanding Fermi gas has been obtained in |9T|.) Thus, the free 
(ballistic) expansion amplifies the scale of correlations. 

It has been shown in [72] that in the presence of a random potential, i.e. 
when the expansion is diffusive instead of ballistic, the picture is different: the 
size of the atomic cloud grows as \/i whereas the (short-range) correlations do 
not get amplified at all. The authors studied the density-density correlation 
function 



C{r,r',t) = {h{r,t)h{r' ,t)) — n{r,t)n{r' ,t) — 5{r — r')n{r,t) , (82) 

where the last term describes trivial correlations, which exist already in a 
classical ideal gas and which are commonly subtracted, in order to isolate 
the non-trivial correlations [90]. Here we only present an approximate ana- 
lytic expression for the normalized correlation function in a "id gas (at zero 
temperature): 



C(r,r',t) 9 _^^n [sin(A;i7'Ar) - {kpAr) cosikpAr)] 



2 



n{r,t)n{r',t) Qs {kpArf ' 

where is the spin degeneracy factor, Ip is the mean free path at the Fermi 
energy, and Ar = \r — r'\. The corresponding expression for the 2d case can 
be found in [72j . 

It is quite remarkable that, while the atomic cloud keeps expanding, the 



normalized correlation function, (83), does not depend on time. The charac- 
teristic length of oscillations remains the same as for the gas in the trap, in 
sharp contrast with the ballistic case, where the spatial oscillation period is 
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growing linearly with time. The peculiar behavior of correlations in a diffu- 
sively expanding gas is related to phase randomization of the wave function 
of a diffusing particle. Eq. (83) implies that, as the expansion proceeds, the 
Fermi gas becomes less "rigid" |72j, in the sense that the relative particle 
number fluctuation increases and approaches the Poisson limit. 



6 Free expansion from a disordered trap 

So far we have discussed two types of problems: The disordered BEC ground- 
state, with the ensuing insulator-superfiuid transition (Sec. 4), and transport 
of cold atoms released from a trap into a random potential (Sec. 5). In the 
present section we consider a different set-up: The trapping potential Vtr{r) 
and the random potential V^(r) coexist in the same spatial region, and the 
atoms are allowed to reach equilibrium in the combined potential. Then, at 
if: = 0, both potentials are switched off and the atomic system undergoes a 
free expansion. We discuss in some detail the case of a strongly anisotropic 
BEC, in a wave guide geometry, and towards the end briefly mention some 
other cases. Our discussion closely follows Ref. |92) . 

The BEC is initially strongly confined in the radial direction, p, by a 
harmonic trap of frequency u±. Weak confinement in the axial direction (z) 
is not essential for our considerations and will be neglected. Thus, prior to its 
release the condensate is in the ground state (we assume zero temperature), 
corresponding to the radial confinement potential |wip^ and a z-dependent 
potential V{z). The latter can be random or deterministic, with character- 
istic amplitude Vq and scale of variation Rq = l/ko- We assume a weak, 
smoothly varying potential, in the sense that Vq << /i and koa± << 1, 
where a± = ^Jliilmuj\_ is the radius of the BEC in the trap. The chemical 
potential /i is assumed to be much larger than huj^. Under these conditions 
the Thomas-Fermi approximation is justified and the ground state density is 

no(p, z) = -^(^- Vi^z) - \muj\p^ . (84) 

At time t = all potentials are switched off and the condensate expands 



according to the equations (44), (45). These equations are to be solved with 
the initial conditions (84) for the density and v = for the velocity field. 
The condensate rapidly expands in the radial direction but, due to the initial 
density modulation, it also develops an axial velocity component Vz{z,t) and 
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the related phase. Upon completion of the first stage of the expansion, at 
time of the order of to = V'^-L? ^^e phase imprinted on the condensate is |93] 



e{z) 



TX 



2huj 



-V{z). 



(85) 



During the second stage of the expansion, for times t » t^, this phase 
imprint can lead to large density mudulation, of the order or larger than 
including the possible formation of matter wave caustics. Since we are inter- 
ested in large effects of this kind, we can neglect from now on the weak initial 
density modulation (its crucial part was to produce the phase imprint 0{z)). 
The second stage of the expansion amounts to linear evolution of the BEC 
wavefunction with the imprinted phase. The wavefunction can be factorized 
into radial and axial parts, \l/(p, 2,t) = ^{p,t)ilj{z,t), and we are interested 
in \ijj{z,t)\^, which gives the density at point z and time t, normalized by 
the radial factor |$(p, t)p. The function %p{z,t) satisfies the linear, time- 
dependent Schrodinger equation, with the initial condition '?/;(z,to) = e'^*^^-*, 
whose solution is 



ij{z,t) 



m f , 

— / az exp 

2'Km J 



im 



J\2 



ie{z') 



it»to) (86) 



where the impressed phase 6{z) can be an arbitrary function of z. The typical 
variation, 6q = Vo/huj±, of the imprinted phase is the important parameter. 
In the case << 1? implicitly assumed in [93j, relative density variations 
remain small for all times, as is clear from expanding the factor expi6'(z') 



in the integrand of Eq (86). Large effects, however, occur in the opposite 



case, 1 << ^0 << (the latter inequality stems from the requirement 

Vq « fi). 

We shall concentrate on the large 6o case, and the purpose is to find the 
form of the relative density \ip{z,t)\'^ at a given instant t. This gives the 
2;-dependence of the density of an expanding atomic cloud, measured with 
a probe beam perpendicular to the condensate axis. The problem is analo- 
gous to the phase screen model in optics, studied extensively by Berry |94) . 
In that model a monochromatic plane wave encounters a thin screen. The 
screen impresses on the wave a phase which may be deterministic or random. 
For strong variation of this phase the wave passing through the screen will 
develop large intensity variations. Observation of the wave intensity at a 
point sufficiently far beyond the screen will reveal a pattern of bright lines ( 
caustics). Although the arrangement we consider differs substantially from 



48 



that in optics (the corresponding matter wave is not at all monochromatic, 
the mechanism of impression of the phase is different, and time assumes the 
role of the spatial axis of propagation in optics), the mathematical treatment 
of the two problems is essentially the same. 

Before discussing the random case, it is instructive to consider a sinusoidal 
modulation 6{z) = Oocoskoz . We introduce the dimensionless variables 
z = koz and t = t/t*, with t* = m/hkQOo, and rewrite (86) as 



^{z, i) = ^ I dC exp [i^o¥'(C, 5, i)] , (87) 

with _ ^ 

(^(C,z,t) = ^^^i^ + cosC. (88) 

The relative density \ip{z^t)\^, initially unity, acquires spatial variations with 
the passage of time. For t <^ 1, relative density modulations are small. 
They grow linearly with t and are oscillatory in z with period 27r. Growth 
of density maxima with time culminates, for 9q 1, in the formation of 
caustics at times t > 1. Caustics are determined by the vanishing of the first 
and second derivatives of the phase y9((^, z,t) with respect to C,. The first 
requirement defines rays of atoms, while the second identifies the points at 
which the emerging rays are focused. The density at these points is found 



by computing the integral in Eq. (87), using the stationary phase method. 



Since the first two derivatives of Lp{C,z^t) vanish, the integral is controlled 

~ 1 /3 

by the third derivative, d'^(p{(, z,t), and has a value proportional to 6q , 

~ 1 /3 

resulting in a large relative density, \'i/j{z,t)\'^ ~ 6^ , with the characteristic 



for caustics aperiodic oscillations [91]. A numerical example is given in Fig. 
5, where the relative density is plotted as a function of z for three different 
values of t (in this example 6q = 30). Caustics are located at 5 = for t = 1 
, and at 5 ?a ±0.7 for i = 2. 
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Figure 5: Relative density as a function of (dimensionless) position z, for 
^0 = 30 and t = 0.5 (top), 1 (middle) and 2 (bottom). (Reprinted, with 
permission, from Phys. Rev. A 80 013603 (2009)). 



Eq. (86) for i/j{z,t) applies also for a disordered potential and the mech- 
anism for caustic formation is qualitatively the same as for the sinusoidal 
potential discussed above. Since the impressed phase 6{z') is now a random 
function, the formation of caustics is a matter of probability. The typical time 
for caustic formation is t* = mRQ/hOo. Using the relations ^mw^a^ = /i and 
^0 = 7rVo/2/ia;_L, the formation time can be written in terms of the experi- 
mentally controlled parameters as 



Quantitative analytic results for the second moment of the relative density, 
1-0(^,^)1"^ = S{t) + 1, can be found in |92) . The calculation is along the 
same lines as for the random screen problem in optics [SS] , where 5* is called 
the scintillation index and it is a measure of spatial intensity fluctuations 
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in the speckle pattern. For a uniform intensity S is clearly zero, and for 
the standard, fully developed speckle pattern 5* is unity. These two limits 
correspond, respectively, to short ( t << 1) and long ( t » 1) times. In 
the more interesting intermediate time regime, when t > 1, 5 is proportional 
to In^o >> 1; signalling the appearance at this time of caustics and the 
associated large density fluctuations. 

Large density variations in an expanding condensate, released from a 
disordered trap, were observed in |17) . In this experiment n/hiJ± = 5.6, 
Rq = 15/im, a± ~ 10/im and the largest density variations were observed 
for Vq = 0.5/i, in a time of flight image at t^oF = 8ms, which is significantly 
larger than l/u)±, see Fig. 6(h). 
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Figure 6: (c)-(j) In situ column density profiles on left and corresponding 
ToF profiles on right: (c),(d) Vq = 0; (e),(f) Vo = 0.3/i; (g),(h) Vo = 0.5//; 
and (i),(j) Vq = l.Ofi. Solid red lines are fits to Thomas-Fermi distributions. 
(Reprinted, with permission, from Phys. Rev. A 77, 033632 (2008)). 

The value Vq = 0.5fi corresponds to a phase amplitude = 4.4 which 
gives the caustic formation time t* = 5ms, comparable to the value of tToF- 
While the value 6q = 4.4 does not lie deep within the large 6*0 regime ( 
the experiment was not designed for caustic observation), the large density 
variations in Fig. 6(h) can most likely be attributed to caustics. On the other 
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hand, a similar experiment performed in [HS] cannot be interpreted in terms 
of caustics, as discussed in [92]. 

So far our discussion was limited to a strongly anisotropic, quasi-lrf BEC. 
Let us briefly comment on some other cases, referring the reader to [92j for 
details. 

(i) 2d geometry. A quasi-2(i condensate is initially strongly confined in 
the axial z-direction (trap frequency w^), with only a weak confinement in the 
radial r^-direction. In addition, there is a random or a deterministic in-plane 
potential V^(r_L). After all potentials are switched off and the condensate is 
freely expanding for a time ~ l/w^, it acquires a phase imprint 9{y±) ~ 



V{Yj^)/hujz. At later times, in analogy with Eq. (86), the planar part of the 
condensate wave function is given by 



$(r±,t) 



m 
27iiht 



d^r'j_ exp 



im 



+ 16 



(90) 



As in quasi-one dimensional systems, caustics are formed for ^ 1 at values 
of the scaled time t ~ 1, and in this regime Eq. (90) can be evaluated using the 



stationary phase method. Unlike the quasi-lrf case, now caustics form a set 
of lines rather than a set of isolated points. The theory of caustic formation 
resulting from a random phase ^(r^) is analogous to the treatment of the 
phase screen problem, which has been studied extensively for two-dimensional 
systems in the context of optics [HI]. In particular, the morphology of caustic 
lines for the random case is discussed in |96| . 

(a) interacting bosons in Id . Here the transverse confinement is very 
tight: hu± is much larger than the characteristic interaction energy so that, 
with respect to transverse motion, all atoms in the trap reside in the 
ground state Xo{p) of the harmonic oscillator, forming a strictly one- dimensional 
system. The axial motion is controlled by an effective one-dimensional Hamil- 
tonian - the Lieb-Liniger model |97| l98j . At some instance the radial confine- 
ment is switched off and, simultaneously, a short potential pulse is applied 
to the system, creating a prescribed phase imprint 6{z) which can be a de- 
terministic or a random function of z. It has a characteristic amplitude Oq 
and scale of variation Rq. The system is beyond mean field theory and the 
whole many-body wave function has to be considered. Initially, i.e. just after 
switching off the trap and making the phase imprint, the axial part of the 
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many-body function is 

N 

1 



^^(^1, . . . zn; t=0) = exp 



$0(^1,... Ziv) (91) 



where $o(-2i? • • • ^n) is the ground state wavefunction, prior to the action of 
the pulse. In the course of expansion, large density variations can develop, 
including formation of caustics. One must, however, account for the momen- 
tum distribution function n{p) of the system of interacting bosons (in the 
mean field approach n{p) = 2tt6{p)). It is clear that a broad momentum 
distribution will impede focusing of the atomic rays and, thus, formation 
of caustics. The width of the distribution n{p) depends on the strength of 
interactions and, roughly, it is determined by the inverse coherence length 
1/C,c- It is shown in [92j that the necessary condition for caustics is C,cni ^ 1, 
where rii is the one-dimensional density. This condition can be satisfied only 
in the weak interaction limit. In the opposite limit of strong interactions 
(hard core bosons) C,c^ ^/ni and the caustics are suppressed. 

(in) Id Fermi gas. Similarly to bosons, one can consider the dynamics 
of fermions, with an initially impressed phase that varies periodically or 
randomly as a function of position z. For a Id Fermi gas the problem is 
similar to that of the previous item. The particle statistics manifests itself 
only through the momentum distribution of particles before the expansion, 
and the Fermi wavelength Ap is taking the place of the coherence length of 
the Lieb-Liniger gas. Since Ap^i ~ 1, caustics are absent from the Fermi gas, 
for the same reason as in the Bose gas with hard-core interactions. The main 
difference between non-interacting fermions and hard-core bosons stems from 
the Fermi surface discontinuity in the momentum distribution. Within the 
approximations of geometrical optics, this discontinuity leads to sharp peaks 
in the derivative of the density (with respect to position z or time t), rather 
than in the density itself. 



7 Conclusions 

Several topics concerning the behavior of cold atomic gases in the presence 
of disorder have been reviewed. This is a relatively new and rapidly pro- 
gressing subject which poses some questions, not previously encountered in 
condensed matter physics. One such question pertains to Anderson localiza- 
tion in a 3d speckle potential, with its rather unusual anisotropic, long-range 
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correlations. It is not yet clear how these peculiarities affect the standard 
picture of Anderson localization and work in this direction has just begun 
|TOl [TT| Ho]. Even such basic quantity as the average density of states is 
not well known for this kind of a random potential and the problem awaits 
careful theoretical investigation. 

In the experimental setup of [TU| [TT] one observes the time evolution of 
an atomic cloud (BEC or fermions) in the presence of a random potential. 
In the long time limit, part of the cloud diffuses away and the rest remains 
localized. To obtain an accurate value for the localization length and for 
the exponent v one has to measure the long tail of the localized part (or 
the very slowly diffusing part, originating form the energy components close 
to i?c)- So far such measurements are beyond reach. Another possibility, 
mentioned in [llj, would be to prepare a collection of atoms with a narrow 
energy distribution and to study its evolution in a random potential. This 
would allow to "isolate" the transition region from the rest of the spectrum. 

One of the most challenging problems in condensed matter physics is the 
combined effect of interactions and disorder. Since non-interacting bosons 
correspond to to a singular, pathological limit, it is clear that the interac- 
tions must be accounted for from the start. Already the Gross-Pitaevskii 
mean field theory for disordered bosons leads to the interesting and difficult 
problem of the Id and 3c/ non-linear Schrodinger equation with a random 
potential. Two of the recently addressed theoretical problems in this context 
are the time evolution of a "matter wave packet" [7H 1751 176] and the flow 
of a BEC through a disordered region [99j. Things get more involved when 
one tries to go beyond mean field. Even the nature of the ground state, and 
the zero-temperature superfiuid-insulator transition, is only qualitatively un- 
derstood HO] . Finite temperature introduces a further dimension to the 
problem and allows for a number of phase transitions. In particular, it has 
been argued that in Id no direct superfiuid-insulator transition is possi- 
ble at finite temperature and that, under increase of disorder, the system of 
weakly interacting bosons must pass through a normal fluid state. (In 'id the 
sharp fiuid-insulator transition is replaced by a crossover). The combined 
effect of interactions, disorder and temperature ( beyond the scope of this 
review) will surely receive much attention in the future. 
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